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INTRODUCTION 


Studies of water resources have relied upon traditional methods of 
data collection* These types of data collection are cumbersome and expenslv(<i 
and result In, at best, a representative sample of a given system* The 
application of remote sensing to the study of the earth's surface has given 
the researcher an Inexpensive tool for the collection of masses of data. 

These data are not a representative sample, but represent the entire area of 
study. The increasing sophistication of the Instruments make possible the 
routine application of remotely sensed data to the study of the environment. 

Throughout 1978 and 1979 an extensive study was made of Utah Lake, by 
the Eyrlng Research Institute, for the Bureau of Reclamation* At this time 
the National Aeronautics and Space Administration was also collecting data 
using the Heat Capacity Mapping Mission (HGHH) Satellite. A grant was 
obtained from NASA to apply the satellite data to the study of Utah Lake. 

The Heat Capacity Mapping Mission (HCMM) was the first of a planned series of 

* 

Applications Explorer Missions that Involve the placement of small, dedicated 
spacecraft In special orbits to satisfy mlsslon^'unlque data acquisition 
requirements (NASA, 1980). The HCMM supports exploratory scientific 
Investigations to establish the feasibility of utlllr^lng thermal Infrared 
remote sensor-derived temperature measurements of the earth's surface. 

The HCiMM spacecraft was launched on April 26, 1978, Into a nearly 
sun-synchronous 620 km circular orbit. Local times of equator crossing are 
2:00 pm and 2:00 am. This places northern hemisphere, mid-latitude over- 
flight times at approximately 1:30 pm and 2:30 am. The selected orbit has a 
16 day ^repeat cycle, with a shorter five day period of approximate repeats 
over some mid-latitude areas. The HCMM orbit covered every area of the 
earth's surface between the latitudes of 85 degrees north and 85 degrees south 
at least once during the day and once during the night within the 16-day 
Interval. 
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OBJECTIVES 


2 • I Surface Tempegature 

The first objective of this research was to develop an accurate model 
for predicting land and water surface temperatures based on the Heat Capacity 
Mapping Mission (HCMM) sensed thermal Infrared radiation* Although a model 
was supplied with the HCMM Users' Guide, the Project Scientist, John C* Price 
(NASA, I9S0) noted that the pref light calibration was "not well understood". 
Llndenlaub and Davis (1978) point out that "...the relationship between 
radiant energy emitted... and temperature. . .Is almost linear over small 
temperature ranges (for e:?ample, between 20 and 3C °C), thus allowing the use 
of a simple linear calibration function for the conversion of radiant energy 
measurements Into temperature values". This statement prompted a decision to 
Investigate the possibility of appylylng a linear model to the data In hopes 
of obtaining an accurate linear correlation. 

In addition to the development of this linear relationship, an 
evaluation of the HCMM supplied model was made. An appropriate offset was 
determined In order to obtain agreement between the model predicted 
temperature values and surface measured values. NASA (1980) reports that an 
offset of approximately -f 5 °C has been applied by some other researchers. 

In an effort to possibly locate any source of discrepancy between 
HCMM and ground measurements, as well as provide additional Information on the 
model, the effects of various atmospheric parameters were also examined. 

2.2 Water Quality 

Another objective of this research was to <^stabllsh and evaluate the 
relationship among HCMM sensed surface temperature and reflectivity on Utah 
Lake and water quality factors Including algae concentrations, algae types, 
and nutrlont and turbidity concentrations. T,he Identification and development 
of measurable relationships would provide a means of employing satellite data 
to yield water quality Information which otherwise Would be collected by much 
less extensive and more costly methods. 


The approach to the water quality* atudy was to analyte and compare 
statistically the measurements made on Utah take by the HGMM Satellite With 
those made by In situ and laboratory testlng> Previous, studies have been 
conducted using Infrared and near-infrared data to detect algae 
concentrations, wind Influences on algae distribution, algae circulation 
patterns, and algae types (Bukata, et al, 1974, Tania, 1978, Plech, et al, 
1978, Matson and Berg, 1981, and Ggan, 1980)» These previous studies employed 
methods such as spectral analysis of algae specimens, multlspectral aerial 
photographs, or positive film transparencies, and algae data sampling to try 
to evaluate the effectiveness of the remote-sensing techniques in Identifying 
algae • 

This study employed the HGMM discrete pixel intensity values along 
with field and laboratory measurements to evaluate the degree of accuracy with 
which the satellite data could Identify algae concentrations and seasonal 
variations* Observations have been made for years on Utah Lake and studies 
have been carried out to determine the types of algae which grow each season* 

Water quality parameters were compared with the HGMM day and night 
Infrared and visible reflectivity measurements. Gomparlsons are on a one-to- 
one basis to determine what relationships exist and with what degrees of 
correlation* Goncluslons have been made about the usefulness of the HGMM 
Satellite data for both directly and Indirectly monitoring the water quality 
parameters* 

2*3 Correlation Images 

Both the day Infrared and reflectance Imagery give a good deal of 
Information about the lake* However, the images separately were cumbersome to 
use* Therefore, another objective was to develop a composite image. These 
were called correlation Images* 

From the , tudy of the day Infrared and visible data It was 
anticipated that unique patterns would become clear In the correlation Images. 
The combinations of Infrared and reflectance Intensities were evaluated In 
order to possibly Identify different phenomena In Various areas of the lake* 
The number of HGMM Images available made It possible to also look for rough 
seasonal variations In the lake* 
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2 .4 Evaporation 

In the arid west evaporation accounts for a major loss of the water 
supply. Calculation of this evaporation Is done with evaporation pans, wind- 
speed functions, and computer models. The pan does not accurately model 
natural water bodies and so the data must be corrected by the use of 
coefficients. Simple wlndspeed functions use local data of humidity, 
temperatures, and wind to estimate evaporation. Most computer models seek to 
eliminate the need for external data by using latitude, daylight hours, 
altitude, solar radiation, etc., to estimate evaporatloii'. 

An objective of this study was to see If remote sensing offers the 
capability of using masses of accurate and comprehensive data In the 
calculation of evaporation. Would there be an Increase In accuracy and a 
decrease In cost by reducing the extensive data collection programs? HCMM 
data was to be applied to empirical models In order to extend the capability 
of the models and enhance their accuracy. The satellite data were used to 
develop a Dalton evaporation relationship and study the effects of wind on 
temperature and hence evaporation. 

Another objective of the research was to study the effects of algae 
on temperature and evaporation through the use of algae grown In evaporation 
pans. For years the effect of vegetation on evaporation has been studied with 
an emphasis on swamp vegetation varying from reeds to trees. Very little has 
been done to address the subject of algae and evaporation. In a remote 
sensing application the effect of algae on water Is Important to know. 

The HCMM data were compared to other climatological data of 1978 and 
1979. The purpose of these studies was to discover the extent to which HCMM 
data corresponded to hydrological data. Pan evaporation was compared to the 
HCMM day and night temperatures and the day/nlght temperature difference. The 
relationship of wind to evaporation and temperature was Investigated. Also 
analyzed were the relationship of the HCMM temperatures to the minimum and 
maximum pan temperatures and the relationship of temperatures in different 
parts of the lake to each other. 


2.5 Thermal Springs 

Several studies have been conducted on the location of thermal 
springs In Utah Lake. Some of the springs around the periphery of the lake 
have been known for years. The late summer and fall HCMM data have been used 
In an attempt to locate areas within the lake where significant thermal 
sources exist. Both day and night thermal infrared data were studied with 
this objective In mind. 

2.6 Groundwater 

It has been suggested that shallow aquifers can affect soil surface 
temperatures. The final objective of this research was to examine the 
feasibility of using the HCMM sensed surface temperatures to locate areas of 
near surface groundwater. The depth of groundwater was defined as the depth 
at which standing water can be found In a shallow hole. The concept of 
thermal Inertia would suggest that diurnal temperature differences would 
possibly correspond quite accurately to the location of near surface water 
bodies. However, the lack of corresponding day/ night HCMM Images constrained 
this research principally to Investigation of absolute daytime surface 
temperatures. An attempt was made to establish the measurable existence of 
such a relationship, as well as determine Its nature. 






a. LITERATURE SURVEY 

3,1 Temperature Calibration 

The HCMM Users’ Guide (NASA, 1980) indicates that there was an 
Instrument calibration problem at launch and that all distributed data are 
susceptible to errors of the order of several degrees (Celsius) . Evaluation 
of data from five satellite passes In May and June of 1978 prompted a decision 
to offset the pre-launch calibration values by -5.5 ^Celsius In order to force 
agreement between the satellite and surface measurements. This offset was 
applied to all standard processed HCMM data. 

The Initial testing of the calibration of the HCMM thermal sensor was 
done on METSA HI, a freshwater lake surface located near the White Sands 
Missile Range. The water surface testing target was chosen over two nearby 
ground surfaces, because of Its more uniform physical characteristics (Barnes, 
1980). 

Further comparisons were carried out for data acquired In October and 
December of 1978, and February of 1979. The results of these comparisons 
suggested an offset of approximately + 5 ^Celsius (in a direction opposite to 
the already applied correction). The more recent data would suggest that the 
prelaunch calibration should not have been modified. 

According to Barnes (1980) the data taken during October of 1978 
Indicated that satellite observed temperatures were substantially cooler than 
surface truth, giving rise to roughly the same order of temperature 
differences as before the adjustment was carried out. During the period In 
question, the water surface In the reservoir used to test the accuracy of the 
calibrated sensor was lower. The surface temperatures, therefore, may not 
have been completely uniform. The surface truth measurements. In turn, may 
not have been truly representative of the large area covered by the HCMM 
pixel. 

As HCMM data are examined In detail by the global scientific 
community. It Is anticipated that this ambiguity will be resolved through a 
substantial number of satelHte/ground truth comparisons (Barnes, 1980). It 
Is also apparent that monitoring of additional parameters describing 
Instrument performance should be considered during the development of future 


missions of this type. Although there Is no evidence that the radiometer 
performance changed over the period of time of the mission, this possibility 
cannot be ruled out based on the apparent change In performance at the time of 
satellite launch. 

Studies are being made Involving the calibration of the HCMM model to 
specific areas. Preliminary results by various researchers have shown that 
the equation Is very sensitive to atmospheric effects, thereby requiring 
careful calibration to any specific area of study. Variations of several 
degrees have been encountered when comparing ground truth measurements with 
predictions from the uncalibrated equation. 

3.2 Water Quality Applications 

In the water resources engineering field, satellites have provided 
especially Interesting and useful information. Landsat data have been used to 
observe and analyze biomass circulation, trophic states and some water quality 
parameters by Welby, et al, (1977) as well as to predict water quality 
according to water color by Chase and Reed (1973) and to detect algal blooms 
by Strong (1972). Satellite data have also been used by Egan (1980) to 
measure turbidity and amrgy absorption properties of water bodies. Gramms 
and Boyle (1971) have studied the feasibility of using selected wavelengths 
for differentiation between algae types. Their research also notes the 
effects of turbidity on reflectance measurements and an increase In algae 
dlstlngulshablllty when phosphorus Is removed from the system. 

Advanced methods of statistical analysis of Landsat data as presented 
by Shlh and Gervln (1980) In the Okeechobee Lake investigation make satellite 
data even more useful to the engineer. Successful comparisons by the 
International Society of Photogramme try, 1980 Working Group 5 (Schott, et al, 
1980) of aerial and surface measurements of several water quality parameters 
also led to the conclusion that large-scale satellite remote sensing 
operations can be very valuable when their application to these measurement 
techniques can be demonstrated. 

A study done by Anderson and Horne (1975) on remote sensing of water 
quality In water bodies In seml-arld climates considered the reflectance 
caused by gas vacuoles produced by Aphanlzomenon flos-aquae, the variations In 
spectral reflectance of blue-green algae as a function of cell activity, 
disruption and death, and the reflectance Of sediment turbidity. This study. 
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together with additional analysea conducted and documented by Horne and 
Wrlgley (1975), Wrigley, et al, (1975), Wrlgley and Horne (1975), and Wrigley 
(1980) concerning lake eutrophication, detection of algae circulation 

;L 

« patterns, and biomass measurements, all provided knowledgeable background. | 

Several studies have 'been made on the growth of algae In Utah Lake [ 

I 

^ » and the parameters which affect It. These Include studies by Squires, et al, 

(1979) and Whiting, et al, (1978) on the competitive displacement factors and 
effects of nutrients and environment on algae growth. Hanson, et al, (1974), 
explain that ERTS-l data were not very useful In detecting algae growth on ■ 

Utah Lake. Research done by Ikeda and Adachl (1978) on the effects of 
nutrients on phytoplankton growth In Lake Blwa Is also applicable to Utah 
.Lake. 

In the visible region the presence of algae Increases the energy 
reflected In the green wavelengths and decreases those In the blue. The HCMM 
band width was large enough to make correction for this unnecessary. Previous | 

studies have preliminarily determined that the presence of algae caused or | 

occurred with Increased water temperatures (Whiting, et al, 1978 and Squires, 
et al, 1979). Several authors have presented detailed lists of Utah Lake 
algae, the cycles of competitive displacement, and the environmental | 

Interactions of the algae (Whiting, et al, 1978 and Squires, et al, 1979). | 

! 

3.3 Evaporation Studies i 

Evaporation Is of such importance in the analysis of water resources t 

that much research has been directed to this effort over the past decades. As 
the amount of data available has been historically small, great effort has 
been made to find suitable relationships that will minimize the need for 
extensive climatological data. One of the first equations developed was the 
so called Dalton Equation which involves an empirical windspeed function. 

Variations of this equation have been developed (Rohwer, 1931) and 
. many forms of the windspeed function relationship have been used In the 

derivation of yet other equations (Penman, 1948, Kohler, et al, 1955, Kohler, | 

^ et al, 1967, and Morton, 1975). The unpublished notes of Dr. A. Woodruff | 

Miller, Brigham Young University, contain many examples. I 

} ' . 
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The Dalton type evaporation relationship was chosen for study because 
of Its applicability to pan data and Its simplicity* This relationship states 
that evaporation Is a product of the wlndspeed function and the vapor pressure 
difference of the water surface and overlying air* 

From published literature, there are several sources of error 
Identified* On small bodies of water a change In wlndspeed has a significant 
effect upon evaporation (Morton, 1975, 1979, & 1980, Penman, 1948)* The 
effect of wind on evaporation decreases with the Increasing size of the lake 
(Kohler, 1955). It Is expected that Utah lake Is of sufficient size to be 
free of the wind effects. Previous studies have given estimates of the pan 
coefficients. 

Vapor pressure measurements can be affected by the vapor pressure 
blanket over the lake* The location at which vapor pressure measurements are 
taken with respect to that of the wind can Induce a "size effect" Into the 
equation* The measurement of vapor pressure should be upwind of the lake In 
air unaffected by the vapor- blanket of the lake. The relationship between 
temperature and saturation vapor pressure Is curvilinear and therefore the use 
of average dally temperatures and dewpoints bias the equation (Kohler, 1955 
and 1967). 

The evaporation model which was to be modified for remote sensed data 
was that of F* I* Morton (1975, 1979, and 1980). This model had been used In 
previous studies of Utah Lake and had been compared to other evaporation 
estimates by Miller, et al, (1980). 

The growth of algae is such a common occurrence In bodies of Water 
that it was thought that a search of available references would be 
particularly rewarding. Very few references on the effect of algae on 
evaporation were found* Rohwer (1931, 1933) suggests that the color of the 
algae Is expected to have a noticeable effect on the amount of evaporation 
because of the change In color of the water* 

There have been studies concerning transpiration and vegetation 
covered bodies (Idso, 1981 and Elsenlohr, 1966). The effect of swamps on 
evaporation has been covered In other studies (Rljks, 1969, Llnacre, et al, 
1970, Llnacre, 1976, and Munro, 1979). 

Size and location are Important factors in describing the evaporation 
of a swamp* Evaporation from the perimeter of an extensive swamp or from a 
long narrow swamp Is enhanced by the prevailing wind* The low humidity and 


the additional energy from the warm winds will Increase evaporation. It has 
been concluded by Rljks (1969) that the amount of evaporation In a swamp 
Increases with an Increase In wlndapeed although the relationship Is 
complicated by canopy roughness. This Is In accordance with the premise of 
the wlndspeed function. 

3.4 Groundwater Studies 

A review of available literature on the use of satellite data for 
determining groundwater depths through correlation with surface temperatures 
yielded very little. Much of the research dealing with HCMH data Is currently 
being done and Is therefore not yet published. The ability of remote sensing 
to detect groundwater, however, has been recognized for more than three 
decades (Reeves, 1975). On conventional aerial photographs, springs and seeps 
are located Indirectly by Identification of particular vegetation types, 
especially In arid regions. 

In expanding the search to Include related material, It was found 
that Rosema and Bijleveld (1978) had done an Investigation using HCMM data to 
estimate the surface relative humidity of bare and scarcely vegetated soils. 
Their procedure Involved estimating the moisture content of the top layers of 
soil. Their conclusion that vegetative cover and soil mineral composition 
both have a significant effect on HCMH values was considered relevant to this 
study. 

Whiting (1976) noted that there Is a direct relationship between the 
measured temperature of bare soil, as well as vegetated areas, and the amount 
of moisture in that soil. Furthermore, although airborne Infrared scanners 
had frequently shot»m variations In soil moisture from a qualitative point of 
view, very little virork of a quantitative nature had been documented. One of 
the primary reasons given for this was the lack of quantitative ground truth 
measurements for correlation with Airborne data. 

Le Shack, et al, (1975), using measurements from a number of 
agricultural fields, obtained a direct linear relationship between soil 
moisture (sampled 10 cm. below the ground surface) and the maximum dally 
surface temperature (Ts max). A sign discrepancy (compared to theory and 
ocher experimental data) was noted, however, and the overall range of surface 
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temperatures examined was only 3 Celsius. Previously, Barthollc, et al, 

(1972) had claimed that Ts max of bare fields decreased as near surface soil 
moisture content Increased. 

Idso and Schmugge (1975) reported that remote sensing of surface soil 
temperature may provide a practical means of assessing soil water status In 
the uppermost few centimeters of bare soil. Also noted was the need to' 
establish the feasibility of obtaining the required temperature measurements 
from satellite remote sensors and more precisely define the depth to which 
water contents could be Inferred from surface temperature measurements. 

Tunhelra and Beutler (1980) reported an attempt to use thermal- 
infrared Images to locate springs and wells through the location of near 
surface groundwater In South Dakota. Two principle conclusions from this 
study were': (1) Bare soil plots with different soil moisture differ In 

surface temperature during the diurnal cycle, and (2) A thick crop canopy 
destroys the apparent surface difference during nighttime hours. They 
Indicated Chat due to a large number of variables (atmosphere, wind, humidity, 
etc.) the differences In surface measurements show much better correlation 
that absolute values. This was attributed to a partial cancellation effect 
when dealing with differences. Also, Tunhelm and Beutler (1980) noted that 
the predawn data was most useful for Identifying shallow aquifers. The need 
for a model relating subsurface soil moisture to surface temperature was 
emphasized In their report. 

Another study by Bonn (1977) gave evidence that grass covered areas 
gave smaller temperature variations than those for bare ground. Ground truth 
measurements showed better correlation to remotely sensed data in vegetated 
areas if taken at approximately two-thirds of the canopy height above the 
ground surface. This would seem to indicate that ambient air temperatures 
have a significant effect on remotely sensed data in vegetated areas and need 
to be taken into account. 

Myers (1970) examined the ability of thermal infrared data to detect 
the presence of shallow acqulfers in an area covered by glacial drift. 

Thermal infrared data were obtained just before dawn during the month of 
August. It was felt that at this time of the year shallow aquifers would be 
cooler than their surroundings and the presence of cool areas in the Imagery 
was expected to indicate the presence of shallow aquifers. The results of the 
statistical analysis Indicated a significant correlation between the thermal 


data and aquifer thickness, thermal dlffuslvlty coefficient, and depth to 
bedrock. Apparently, the correlation with depth to water and depth to the top 
of the aquifer were not statistically significant. 

: \ 
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4. 


STUDY AREA 


4.1 Utah Lake 

This HCMM hydrological study was carried out using Utah Lake (Figure 
4.1) as the primary target. Utah Lake Is located In central Utah, 
approximately 40 miles south of Salt Lake City and la the primary freshwater 
lake In the Jordan River drainage basin. The lake Is one of the largest 
freshwater lakes west of the Mississippi River, covering approximately 385 km 
(95,000 acres) at "compromise level" (1368.35 m, 4489.34 ft.). The compromise 
elevation Is the established level to which the lake Is permitted to rise 
before free flow through the outlet works Is allowed. 

Although the lake covers over 25% of the valley floor. It has an 
average depth of only 2.4 m (Fuhrlman, et al, 1974). Utah Lake Is a shallow, 
unstratified- eutrophlc, arid region lake with approximately 50 surface 
Inflows and one surface outflow. The water Is highly turbid with Secchl depth 
readings averaging 24 cm and ranging from 12 cm to 50 cm. The lake basin 
receives Inflow from numerous mineral springs within and around the periphery 
of the lake and as a result the water has a high carbonate and sulfate 
content . 

An accelerated eutrophication process has occurred In Utah Lake. 

When the Mormon pioneers first settled in this area the Utah Lake fish 
population was dominated by a species of cutthroat trout (Salmo clarkl) which 
was adapted to the already slightly eutrophlc conditions of the lake. Due to 
the use of Utah Lake and tributaries for Irrigation and waste disposal and the 
Introduction of other fish species Into the lake, the trout have become 
extinct (Whiting, et al, 1978). 

The growth of algae in Utah Lake has Increased with the Increase In 
nutrient loads caused by sewage disposal Into the lake and the return flows 
from Irrigated fields, pastures and corrals In the drainage basin. Algal 
blooms on the lake during late summer and early fall are common. The large 
algae concentrations In Utah Lake cloud and discolor the water and produce 
obnoxious odors and decomposing mats of organic matter on the surface and 
around the shoreline. Perhaps the most serious concern Is the depletion of 
dissolved oxygen In the lake caused by the decaying algae. All of these 
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factors combine to reduce the attractiveness of the x^aters of Utah Lake for 
use as a domestic or industrial water supply and as a site for fishing, 
swimming or skiing. 

The lake is considered an ideal size for the satellite Imagery study, 
covering approximately AO X 70 HCMM pixels. A large water body has several 
distinct advantages for a remote sensing hydrology study. Several of the more 
Important advantages are homogenlty, thermal stability, minimal relief, and 
size. Jackson (1980), in attempting to use HCMM imagery. Indicated that it is 
extremely difficult to locate a nine pixel sized field. 

Utah Lake is a vital part of the Central Utah Project being developed 
to fulfill the need for water in these valleys. Measurements of water quality 
parameters in Utah Lake have been made for several years, including 1978 and 
1979. The existence of the large number of recorded measurements made Utah 
Lake the ideal location for analysis and comparisons of ground-truth 
measurements with the data returned by the HCMM Satellite. 

A .2 Southeast Utah Valley 

The area immediately to the southeast of Utah Lake was chosen to 
study the correspondence of HCMM intensities to near surface groundwater. 

This area is chiefly agricultural, but is broken by large tracts of land too 
poorly drained to be effectively farmed. These poorly drained, "high 
groundwater” areas are ideal for the groundwater portion of this study. 
Numerous areas cover several pixels specially and due to their close proximity 
to Utah Lake, are relatively easy to locate. The actual area examined extends 
from the foot of West mountain, in an easternly direction, to the residential 
area of Sprlngvllle, Payson and Spanish Fork. The southern boundary is 
roughly between Payson and the southernmost extremity of West mountain, while 
Utah Lake marks the northern limit. The study area is generally quite flat, 
largely between 5 and 20 feet above the surface of the lake. 
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HCHM DATA 


5 .1 Instrumentation 

The HCMM satellite carries a single sensor ^ the Heat Capacity Mapping 
Radiometer (HCMR). It acquires high spaclal resolution daytime reflected and 
thermal Infrared data. The HCMR has a geometric Instantaneous field of view 
(IFOV) of 0.83 mllllradlans, high radiometric accuracy, and a wide coverage 
swath on the ground. The Instrument Is a two channel scanning radiometer, one 
spectal channel covers the reflectance band from 0.5 to 1.1 micrometers, while 
the other channel views the thermal Infrared band between 10.5 and 12.5 
micrometers. The two channels thus provide measurements of reflected solar 
and emitted thermal ener$;y, respectively. The overall measurement accuracy In 
the thermal channel Is approximately 0.5 *^C at 280 °(C. 

The spaclal resolution of the infrared channel Is approximately 600 
by 600 meters at nadir. These values are masked by data processing, which 
generates registered data at a 481.5 meter pixel size. The swath of data 
coverage along the track Is approximately 716 km wide. The scale of Imagery 
Is nominally 1: 4,000,000. The HCMM satellite officially terminated Its 
mission In October of 1980, having provided reliable data chiefly during 1978 
and 1979. 


5.2 Data Products 

For convenience master output tables were used which consisted of 
digitized raw counts from both channels 1 and 2 of 8~blt (0-255) Indices. 

This yielded the result that a count value from any of the HCMR data tapes 
always represented the same value of radiance or temperature. These data were 
used to produce geometrically corrected Images of day and night temperature 
and day radiance Intensities. 

HCMM Satellite data employed in this study were supplied by NASA as 
positive and negative prints of thermal and reflective Intensities, 
represented by grey tones, together with computer compatible tapes. A total 
of 23 HCMM daytime Imajjss were available for the study area. Five HCMM night 
Images were also available. Table 5.1 lists the dates of Images used. 
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Utah Lake Location 
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The computer tapes were used to generate line printer maps on which 
Utah Lake was located In order to match particular pixels with measurements 
made on the lake Itself at about the same time. The HCMM scene Is about 700 
km square. From this It was necessary to locate and extract the Utah Lake 
window. The starting row and column of this window within the full Image were 
established. The distance was measured to the corner of the window on the 
NASA print. This was then used to calculate the number of rows and columns. 
The VAX computer system was used to locate and extract the Utah Lake windows. 

HCMM digital thermal and reflective data were printed on a line 
printer as a two dimensional array In a Fortran 13 format. Locating the Utah 
Lake study area was rarely difficult, due to the characteristic difference 
between land surface and water surface radiation values (Figure 5.1). 
Differences were on the order of twenty to fifty Intensity units (7 to 16 °C) 

Sampling Stations Location 

Location of corresponding township and range lines was necessary 
because most available surface data was referenced to the State Township and 
Range System. While the lake body was easily located, the scale was severely 
distorted due to the digitized format. It appeared as roughly a 3 horizontal 
(east-west) to 1 vertical (north-south) proportion as compared to actual 
ground distances. The daytime lake map was skewed approximately II degrees, 
apparently as a function of the flight path of the HCMM Satellite. The 
nighttime images exhibited a skew of approximately 15 In the opposite 
direction. 

Township and range lines were located by using geometric 
trlangulatlon and scaling ratios, employing several prominent features located 
around the lake as reference points (Figure 4.1). The east-west protrusions 
were much more prominent on the HCMM digital Images than were north-south 
protrusions. Due to this amplified prominence, east-west protrusions were 
used wherever possible In the trlangulatlon procedure. 

Having established the most prominent trlangulatlon points, 
alternating section lines, both north-south and east-west, were scaled 
vertically and horizontally across the study area. The sealing was done by 
using U.S. Geological Survey quadrangle maps. These section lines were then 
transferred to an overlay of the Utah Lake line printer maps, using the 
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Figure 5.1. HCMM Sensed Land Surface Emissivity vs. Water Surface 
Emissivity (Goshen Bay Area) 


19 


distance ratios as determined from the quadrangle maps. The non~scaled 
section lines were then drawn exactly' between the scaled lines. This raethod 
proved to yield sufficiently accurate overlays of the study area. 

Overlays of Utah Lake were prepared for each format and the sampling 
points identified on them. These overlays were then placed on the HGMM data 
printouts to identify the pixel Intensity values which corresponded with the 
local sampling stations for any particular day. Tables containing the pixel 
Intensities together with the other measurements made at the sampling sites I 

for each day when data were available were uiade. 
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GROUND TRUTH DATA 


6tl Sampling Stations 

During the period when the HGMM was collecting data> several studies of 
Utah Lake were also carried out* These studies, sponsored by the U.S. Bureau 
of Reclamation, produced the data with which the HCMM Information was 
evaluated. The Water Quality, Hydrology and Aquatic Biology (WHAB) study was 
conducted by the Eyrlng Research Institute, of Provo, Utah and consisted of an 
extensive sampling and evaluation of Utah Lake waters at 15 stations 
throughout the lake. The locations of these sampling stations are indicated 
on Figure 6.1. The prefixes Indicate the area of the sample stations, l.e., 
UL'-Maln body of Utah Lake, GB-Goshen Bay at the southern end of Utah Lake and 
PB-Provo Bay on the eastern side of the lake. 

Transect sampling was done on September 28 and November 21, 1978 and 
April 25 and 27 and July 16 and 17 of 1979. These transects were established 
In separate parts of the lake In order to determine whether the biota at a 
point In time were homogenous throughout the lake or shoved local variations. 
Samples were taken quarterly so that seasonal variations In the biota could 
also be determined. Figure 6.2 shows the locations of these transects and 
sampling points. Information on the WHAB studies included herein Is from 
Rushforth, et al, (1980 and 1981). 

6 .2 Temperature Measurements 

Water temperature, being a relatively important water quality parameter, 
was routinely Included in the water quality data. Examination of this data 
revealed that fourteen stations had data pertinent to this study. These 
stations are chiefly In the central parts of the^lake, thus reducing the 
possibility of the shoreline Influencing the radiation measured by the HCMM 
sensor. Data from the stations that were In closer proximity to the shore, 
primarily the Provo Bay and southernmost Goshen Bay locations, were analyzed 
carefully to insure that no land surface overlapped the pixel boundary. WHAB 
data are tabulated In Merritt, et al, (1981). 
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The other Important source of temperature calibration data was a study 
done by the Utah State Division of Wildlife Resources. This project Involved 
the capture of fish samples for research dealing with the aquatic biology of 
the lake ecosystem. 

Water temperature and location were routinely tabulated , along with 
other data, for each sample taken. The majority of this data were taken from 
a small b<u,it, but due to the objectives of the project, were generally quite 
close to shore. Because of this, extreme care was taken to minimize the 
probability of Including land surface effects In the HCMtJ values. In cases 
where It was felt that significant land surface war/ Included In the pixel, the 
Intensity value from the adjoining (lake side) pixel was used. In cases where 
the location appeared to fall between pixels, the arithmetic mean of the 
surrounding pixels was used. 

While It was felt that the surface measurements were of relatively 
similar accuracy, one other factor needed to be taken Into account. The HCMM 
Satellite was designed to pass overhead at the warmest part of the day, 
approximately 1:30 p.m. The surface data, on the other hand, was generally 
collected throughout the course of the day. For this reason. It was 
determined that a certain degree of fluctuation should be expected In the 
data, simply to account for the daytime surface temperature fluctuation of the 
study area. 

Examination of all available day/night HCMM images with a 12 hour 
Interval of separation showed that characteristic lake areas had roughly a 3 
(9.2 HCMM intensity units) diurnal temperature variation. It was also 
determined that this value was most likely slightly low. This determination 
was based on two facts. First, the nighttime temperature was sensed by HCMM 
at approximately 2;30 a.m. The coolest portion of the night, on the other 
hand, occurs slightly before dawn. Therefore, the nighttime HCMM reading was 
felt to be higher than the true "nighttime low" temperature. Second, all 
available day/night corresponding Images were during the months of August and 
September. Since calibration data was used for periods from March through 
November, It was felt that the seasonal fluctuation in diurnal differences was 
not truly represented. The August-September period in fact, had the smallest 
diurnal fluctuation of the study period. Based upon the above observations 
and calculations. It was determined therefore, that the measured data would 
likely have a fluctuation of around 2 °C. 
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6«3 Alj^ae Measurements 

Phytoplankton samples were collected in one of two ways. If the water 
contained large numbers of algae, lake water was collected directly by 
Immetl'slng a 15 liter bucket beneath the surface of the water and lifting It to 
the surface. A subsample was then removed from the bucket and placed In a 
Nalgene container. These containers were returned directly to the laboratory 
for analysis. Samples collected In this manner were designated total plankton 
samples. 

When the algal concentration In the water column was too low for direct 
water analysis, known quantities of water were passed through a plankton net 
and the algae were collected In attached 30 ml vials. These samples were 
designated as net plankton samples. When these samples were collected 
corresponding nannoplankton samples were usually taken by catching a portion 
of the water that was passed through the net In a Nalgene Jar. These were 
returned to the laboratory for analysis. 

Standing crop measurements were made by counting algae using Palmer Cell 
counting chambers. These chambers were placed on the microscope and the algae 
present In the field of vision were counted under 400x magnification. The 
slide wa» then moved to the adjacent field and the algae In that field 
counted. In this manner the slide was moved from end to end and the number of 
algae counted In this transect was recorded. The number of organisms present 
In the water of the lake was then calculated by multiplying the counted 
quantities by determined factors. 

Several statistical analyses of the data were performed, Including, 
analysis of diversity using the method of Shannon and Weaver (1963). 
Calculation of Importance values (Ross and Rushforth, 1980) was also performed 
as well as several other tests. 

6.4 Climatological Measurements 

Climatological data needed were readily available. Humidity data was 
collected and tabulated during the study by the Bureau of Reclamation. 

Humidity data was taken on the east side of the lake. This is on the downwind 
side with respect to the prevailing winds and as such may Induce some bias. 

Evaporation Is measured at two sites along the lake, the Bureau of 
Reclamation (USBR) pan at the Provo Airport and the Utah Lake Lehl weather 
station. Both sites have wind data and Lehl has pan and air temperature data. 
Both pans are close to the edge of the lake. Wind data from Provo Is at both 


a height of 11 tS m and at pan height* tehl wind Is at pan height. The 
velocity of wind over water changes relative to that on land and Increases 
depending upon the fetch. As both sites are on the downwind side of the lake, 
wind velocity measurements correspond closely to that of wind over water. 

Evaluation of the HCEMH temperature and evaporation data yielded fourteen 
days of HCMM day infrared data and five day/nlght infrared pairs that could be 
used In the analysis. HCMM temperatures were used to calculate the saturation 
vapor pressures of the surface of the lake. Maximum air temperatures from the 
Lehl station were used In calculating the prevailing vapor pressure. 

6.5 Groundwater Measurements 

An extensive search through Federal, State, County and local sources 
revealed little reliable quantitative near-'surface groundwater data. The 
Bureau of Reclamation study Included extensive groundwater sampling, but was 
primarily concerned with deeper aquifers and provided no useable data. 

Similar situations were encountered with the U.S. Geological Survey, the Utah 
State Geological Survey, and the Utah State Engineers office. The only 
quantitative data available was from Utah County Watermaster records of 
percolation testing. County authorities felt that some of the data was quite 
reliable. Part of the records, however, were from tests not supervised by 
county officials and such tests were marked accordingly. The lack of 
reliability of these unsupervised tests caused the requirement for official 
supervision to be Instituted In 1979. Nevertheless, all data, whether from 
supervised or unsupervised tests, were tabulated and used in the development 
of the quantitative groundwater relationship. The possible differences 
between predictions based on supervised and unsupervised tests were 
Investigated. 

After falling to find a reliable and extensive source of depth to 
groundwater data, a program of field Interviewing was Instituted. Local 
residents, farmers and ranchers chiefly, were Interviewed in an effort to 
recreate the groundwater history of the area throughout the study period. 

This method met with limited success for several reasons. Individuals were 
generally able to furnish Information concerning only their own fields, or 
fields that they worked. There was an extremely high number of . 
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contradictions, some with differences of up to several feet, In neighboring 
and similar areas. Feeling that this procedure was not yielding useable data,^ 
It was discontinued. 

Because of the limited data base from which to develop a statistically 
significant quantitative relationship. It was decided to also explore the 
possibilities of proving the existence of such a relationship from a 
qualitative point of view. The search for quantitative data had produced a 
number of sources of data capable of being used In a broader and less precise 
study. The most valuable was a "wet season water table” map developed by Ben 
Everltt (1981), formerly of the Utah State Geological and Mineral Survey. 

This map was based on existing data obtained from reports and files, both 
published and unpublished, as well as Interviews with Various government 
agencies and consulting firms. With the exception of a Bureau of Reclamation 
study, most of the data were one time readings of water levels In temporary 
shallow bore holes or trenches. Meagurements were generally from the 1960's 
and 1970' s. 

The Bureau of Reclamation study consisted of a detailed drainage 
Investigation In Southern Utah Valley and Goshen Valley. This investigation 
was done as part of the planning phase of the Bonneville unit of the Central 
Utah Project. More than one hundred shallow water table observation holes 
were monitored over the fifteen year period from 1953 to 1967. 

Care was taken to use onl 3 i' data from the shallowest boreholes and 
excavations available. This was done In order that the Influence of the 
plezometrlc heads of the deeper, confined aquifers upon the water table could 
be minimized. Toward the center of the valley, artesian pressures at very 
shallow depth have been reported, in two cases within ten feet of the surface 
(Everltt, 1981). 

Another significant source of regional data was the "Soil Suitability 
for Septic Tank" study done by the U.S. Soil Conservation Service (SCS). This 
study provided detailed maps of locations In the county that experienced 
seasonal water table depths of four feet or less. Another map, developed by 
the SCS In conjunction with Utah County, showed poorly drained areas, as well 
as those subject to periodic flooding. 

The above sources, combined with Interviews with many local reslddhts, 
water board members, and county officials having first hand knowlege of 
specific areas, were used to develop overlays showing areas that experienced 


seasonal variation in groundwater depth. Special attention was focused on 
nonhomogeneous location; e.g. , a deep water table area surrounded by shallow 
water table areas, and vice versa. 
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ANALYSIS METHODS 


7*1 Sum of Squares of Residuals 

One of the primary statistical concepts used for the determination of 
the comparative correlation of predicted data to actual ground truth was the 
"sum of aq*Hares of residuals" (SSR). This value Is determined by subtracting 
the predicted temperatures from the measured ground temperatures, squaring the 
Individual differences, and obtaining the sum of these squared terms. The 
model yielding the smallest SSR Is chosen as the model best predicting the 
actual data. This Is a widely accepted and commonly used method of comparing 
the relative "fit” of different models to actual values (Walpole and Myers, 
1978). One important stipulation, however. Is that the two models must each 
contain the same number of predicted values, and these In turn must have a one 
to one correspondence with measured or "true" values. 

The effect of squaring the differences before summation possibly 
warrants further explanation. The terms are squared chiefly for the purpose 
of making all summed values positive. If this were not done, positive and 
negative values would cancel, and the comparative value of the respective sums 
would be lost. 

One other effect of using the SSR as the chief Indicator of fit Is 
the exaggerated significance given to points lying significantly off of the 
predicting line. Since the difference Is squared, a single "outlying" point 
can have a significant Impact on the SSR. For this reason, data used to 
compute the SSR should be examined very closely for possible points that could 
unduely Influence the results. 

Possibly the most Important tool for determining the goodness of fit 
of a specific statistical model Is the analysis of residuals. The residual Is 
the difference between actual and predicted values. Once the "best fit” model 
has been determined through use of the SSR, the residual values are plotted as 
dependent variables. Residuals can be plotted against parameters such as data 
location or time of acquisition. Any pattern emerging from a residual plot 
warrants close examination. 
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Although the SSR value can give a good Indication of the best fitting 
model, It Is not sufficient to determine If one model is statistically better 
than another, unless one model Is a special case of the other. As this was 
not the situation with the models being compared, another test was needed. 

The available data was randomly divided Into two portions. One set of data 
was used to determine the best linear model and then this same data was used 
as the test data base for the NASA model. The other set of data was used to 
determine the calibration coefficients for the NASA model, then used as the 
test data base for the linear model. This rather awkward arrangement was 
determined necessary In order to meet the mathematical requirements for a 
valid statistical test. 

Each model was then applied to Its respective data base, aiid the 
Individual SSR values were determined. Each SSR was then divided by the 
number of data points used In Its generation. The ratio of the two values 
obtained In this manner then met the requirements stipulated to follow a 
statistical F distribution. By forming an appropriate set of hypotheses and 
comparing the value obtained (the ratio) against tabulated values of the F 
distribution, the superiority of one model over another can be either proven 
or dlsproven (from a statistical point of view). Again, It was not the Intent 
of this study to develop a better model, but rather to determine the 
feasibility of developing a computationally simpler model without sacraflclng 
accuracy. 


7 .2 Linear Regression 

The best linear model was detexrmlned through the use of simple linear 
regression. This Is based on determining the line with the smallest SSR. The 
regression analysis, as well as most all data handling and manipulation, was 
done with the aid of Mlnltab II, a statistical software package developed by 
Professor Thomas A. Ryan Jr., of the Statistics Department at Pennsylvania 
State University (Ryan, et al, 1976). 

Mlnltab II Is a general purpose statistical computing system. It Is 
* easy to use, very flexible, and fairly powerful. It la designed primarily for 
small to moderate sized data sets, which can be stored In main memory. 

Mlnltab has been found especially useful In exploring data In the early phases 
of analysis, for plotting, and for regressions. 


7 *3 Goi?g(ilatii.oni> 

tn ogdttg to ovaluHto tho tolatloushipis tUot oxtst botwoon th« varlout 
raoxMUgQmoot8 used In this (itudyi «Hiv«va]. ngoumptlout vtg« w«do* Thoii« tnoludo 
tht g«lAtiv« Indopondonot of «aol\ of tho Mioturomontt » « otmUor d«gr«« of 
ncmigacy of oaoU off th« mtasurowonts and a notwal dlatrlbution of oach of tha 
data gata< 

Tho aisumption of Indapandonoa tinpXlos only that no talatlonahip vas 
aaauaod until It could bo pvovoui fot axawpla, waa water tomparatura dapandaut 
on the algaa concavit ration or was algaa coucantratlon dapandaut on tUa watar 
tomparaturat In ovary caaa tUo corralatlou off the data nets aa an indapandant 
bivariate function f(x,y) was aatabllaUad boffora performanca of any linear 
ragreaslon analyaon which roiiulra the assumption off dapandenca of one variable 
on the othar« the determination of the dagraa of correlation of the data 
sets I together with Information gathered from literature on previous research 
concerning the meaeur omenta In queatloui provided a key to the dependent 
relatlonahlps that existed. 

the assumption of a similar degree of accuracy off the measurements 
allowed the use of one-to-one correlation off the data sets. In the cases 
where this assumption appeared to be unjustified or questlonable» either 
adjustments were made and recorded as part off the evaluation procedure, or no 
correlation was established because of the Incompleteness of some data sets. 
One principle adjustment wlilch was made was the use of the log of the plankton 
counts, rather than the counts themselves. In order to evaluate data sets of 
similar magnitude and to Increase the linearity^ of the sets. 

The assvimptlon of a normal dtatrlhution of the measurements In the 
data sets was necessary dne to the often limited number off common points (less 
than 30) in the sets. This type of assumption la common to statistical 
analysis since It makes possible the substituting for a very difficult 
equation for evalviatlng the degree off correlation with a much simpler equation 
using the t-dlstrlbntlon. Since the t-dlstrlbntion Is similar to the normal 
distribution in tl\at they both are syimetrlc about a mean off aero and ate 
bell'"Shaped, there is no inconsistency In the simplification. Also, since the 
t-distrlbutlon Is more variable, owing to the fact that the t values depend on 
the fluctuation of two quantities (the sample mean and variance), the accuracy 
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of the normal approximation ia increased by the simplifying procedure * As the 
sample size increases to infinity the standard normal distribution and the t- 
distribution approach equality. 

The correlation procedure used in this study assumed that both the X 
and Y sets in each comparison were random variables and the measurements were 
observations from a Joint density function f(x,y). According to Walpole and 
Myers (1978) it is often ossumed that the conditional distribution f(x|y) of 
Y» for fixed values of X» is normal with mean y(Y|x) " 0( d- 3 x and variance 
of^ (Yjx) ■ ax^ and that X is likewise normally distributed with mean u(x) and 
variance (x). 

It was important to understand the physical interpretation of the 
correlation coefficient and the distinction between correlation and 
regression. The straight line given by y(Y|x) ^ a + 3X la the regression 
line. Since r^ - $^0x*/ay^ (Walpole and Myers, 1978), the value of the 
correlation coefficient r is zero when @ >■ 0, which results when there is 
essentially no linear regression; that is, the regression line is horizontal. 

A value of r equal to -fl Implies a perfect linear relationship with a positive 
slope, while a value of r equal to results from a perfect linear 
relationship with a negative slope. 

7*4 Color Graphics 

One of the most valuable methods used in this study to Interpret the 
HCNM Satellite data was the preparation of color representations of the pixel 
intensities on Utah Lake. Those color pictures were made on the Tektronix 
4027 color terminal in the BYU Civil Engineering Graphics Laboratory. The 
IICMM data tapes were read and stored on the DEC~10 computer and an Interface 
was used to read these files into the Tektronix 4027. The range of pixel 
intensities on the lake was divided into seven smaller ranges and each of 
these had a color assigned. All pixel intensities higher or lower than %'hose 
assigned were represented as black. Because of the normal difference between 
the water and land pixel intensities of 20 to SO units the ranges could be 
selected so that the lake appeared in color on a black background. The colors 
were assigned so that they ranged from dark blue to red as the intensity 
measurements Increased from the lower ranges to the higher. 
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These color pictures were especially useful Ip locating the extent 
and relative concentrations of algae In the lake waters and In Identifying 
blooms within the algae fields. They also proved to be very valuable in 
Identifying regions In the lake where some additional thermal energy appears 
to be present. When Insufficient data were available to do reliable 
statistical analyses of the measurements, the color pictures sometimes 
provided the only means to evaluate data. In almost every case where 
statistical analyses could be compared with the color pictures, the pictures 
were found to be truly representative of the concentrations or other 
measurements In question. 

In order to generate the Infrared/reflectance correlation images a 
two dimensional graphics program developed In the BYU Civil Engineering 
Department was modified and used. The amount of color that could be displayed 
was limited and therefore the degree of refinement to be applied to the 
combination of the Images was reduced. The pixel values of the day Infrared 
and visible data were read Into the program which split each data set into 
three ranges; high, middle, and low. The upper and lower limits of these 
ranges for the lake were established by Inspection of the lake Image printout 
In order to exclude edge effects at the shoreline. 

The program compared each pixel Intensity of the Infrared Image to 
the range limits and determined If the pixel value was In the top» middle, or 
bottom third of the lake Image Intensities. Each pixel value for the visible 
Image was also read and assigned to the top, middle, or bottom third of the 
day visual Intensity range. The Infrared and day visible Images were then 
compared as shown In Table 7.1. From this comparison a new Image was formed 
In which the pixel values were the numbers assigned (columns 3, Table 7.1) 
from the comparison. The image was then read by the graphics program, 
displayed on the screen, and a picture was taken of the image. 
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Table 7.1. Numbers and Colors. Assigned In the 
Comparison of Day Infrared and Visible Images 
to Produce the Correlation Pictures. 


Day 

Infrared 

Day 

Visible 

Number 

Assigned 

Color 

Assigned 

H 

H 

7 

red 

M 

L 

6 

pink 

M 

H 

5 

yellow 

M 

M 

4 

green 

H 

L 

3 

It. blue 

L 

H 

2 

blue 

L 

L 

1 

die. blue 

H 

H 

0 

black 

L 

M 

0 

black 
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TEMPERATURE CALIBRATION 


8tl Theory 

The dominant energy flu>;es controlling the earth's surface 
temperature are Insolation^ during daylight hours, and radiative and 
convective cooling, which occur throughoitt the 24 hour period. Solar Input Is 
concentrated In a wavelength Interval centered at 0.5 micrometers, while 
radiative cooling to space occurs In a broad, longer wavelength spectral 
Interval of from 6 to 30 micrometers. Alternation of the two fluxes causes a 
diurnal cycle of surface temperatures. The HCMM Is designed to sense surface 
Infrared radiation values at thq peak of the cycle (1:30 pm local time) and 
near the early morning minimum (2:30 am local time). 

An Important concept In the use of thermal Infrared Imagery Is the 
differentiation between reflectivity and smisslvlty. The range o‘k daily 
surface temperatures Is a function of a physical property of a material to 
absorb and re-emlt the suns energy. Since this energy la conducted to a depth 
of two to five centimeters, the measurement of surface temperature variation 
provides Information on a finite thickness layer of material being observed, 
and the sensed energy Is a function of the type or state of material. This 
contrasts with measurements of reflected energy, which only provide 
Information about the surface of the observed material. 

Interpretation of surface temperature measurements obtained by HCMM 
are complicated by many factors, most of which are associated with the earth's 
atmosphere. Several of the more Important factors are; cloudcover, surface 
winds, evaporation, water vapor, surface emlsslvlty, directional 
characteristics of reflecting materials, the prior temperature history of the 
observed material, topography, and surface variability within the area of a 
HCMM pixel. In a general sense the complete energy and moisture budget at the 
earth's surface must be considered. The effects of some of these parameters 
were examined. 


A spaceborne sensor must be ’able to “see" through the atmosphere to 
the ground without attenuation* Radiation may be lost In Its passage through 
the atmosphere as a result of scattering from and absorbtlon by suspended 
particles* The atmosphere's minor constituents, such as water vapor, carbon 
dioxide, and ozone create the major absorption bands. 

The relative Importance of the different absorbing molecules depends 
on the radiation path. A satellite located sensor receives radiation that has 
had to pass through the entire atmosphere. The HCMM thermal sensor waveband 
of 10.5 to 12.5 m corresponds to a rather pronounced window In the far 
Infrared region. 

8.2 Offset Calculations 

8.2.1 Preliminary Analysis 

The available measured temperature data were tabulated, along with 
corresponding dates, atmospheric parameters and HCMM radiation Intensities. 
Also tabulated was a code identifying the region of the lake from which the 
data had been obtained. The lake was divided Into four general regions, 
respectively (1) North Lake, (2) Middle Lakej, (3) Provo Bay and (4) Goshen 
Bay. Although a large amount of data was tabulated for Provo Bay, much of It 
was during the late summer period. At this time, a large percentage of the 
Bay Is essentially marsh, with significant expanses of water plants covering, 
shadowing or otherwise concealing the water surface. Feeling that the 
vegetative cover could possibly have a detrimental effect on the calibration, 
this area was carefully examined and much questionable data was eliminated 
from the tabulation. The tabulated data were then entered Into Mlnltab for 
manipulation and analysis. 

The first operation performed was to plot the surface measured 
temperatures over the curve predicted by the HCMM model. As can be seen In 
Figure 8.1, the HCMM curve seemed to fall somewhat below the majority of the 
ground truth data. In order to further clarify this, a plot was made of the 
surface measured temperatures against the HCMM predicted temperatures (from 
Users' Guide equation) in Figure 8.2. If the two were in reasonable 
agreement, the data would then be expected to fall generally on a 45*^ line, 
starting at the origin and bisecting the two axes. As can be seen from Figure 
8.2, this was not the case. The data fell generally to the left of the line. 
Indicating that the HCMM equation was predicting the surface temperatures to 
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be lower than they actually were« This appiiars to agree with the conclusion 
of Jackson (1980) that temperature data obtained from the HCMM were 
consistently lower than temperatures measured with a low altitude thermal 
scanner. 

Offsets were then added to the HCMM predicted values and the SSR was 
calculated for each of the applied offsets. Figure 8.3 is a plot of the SSR 
value verses the applied offset. As can be seen, the required offset to yield 
the smallest SSR, and thus the best fit, was about 4.6 °C. Figure 8.4 shows a 
plot of ground measured temperature against the HGMM predicted temperature 
with the 4.6 °G offset, along with a + 2 °G envelope. As can be seen, a 
significant number of the data points lie outside the envelope. 

In an effort to determine the possible causes of these deviations, 
several plots were made. A plot of the residuals verses lake location yielded 
no noticeable pattern. The plots of residual verses ground temperatures 
(Figure 8.5) and residuals verses months of the year (Figure 8.6) did yield 
usable results. Examination of Figure 8.5 showed that a group of data in the 
10 to 20 °G temperature range had negative residuals. Examination of Figure 
8.6 showed that data collected in the month of May, (5 month) also had 
negative residuals. However, plots of the residuals versus location, month, 
and ground temperature (Figures 8.7, 8.8 and 8.9, respectively) for only 
"same-day" data ( with a presumed 5 offset) did not show the deviant 
points. This indicated that the data in question was collected on a day other 
than the day of the HGMM overflight. 

In order to further Isolate the data in question, the WHAB data was 
separated from the Fish and Wildlife data, and these sets further subdivided 
into delta day categories. Delta day categories are groups of ground surface 
data taken a specified number of days apart from the HGMM overflight. Figure 
8.10 indicates that two groups of the WHAB data fall noticeably below the 
expected 45° line (offset of + 5 °C added to HGMM temperatures). One group 
corresponds to a delta day of four days, the other group to a one day 
difference. 

Upon detailed examination of the data, it was found that the Fish and 
Wildlife data collected during May 1978 were all in the 10 to 15 °G range, and 
all had large negative residuals. The first three plots, therefore, were all 
Implicating the same data. An investigation of the time Intervals further 
substantiated these Implications. The HGMM overflight corresponding to this 
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Figure 8.3. SSR vs.- Applied Offset: All Tabulated Data 
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Figure 8.6. Residual vs. Month: All Data 
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Figure 8,10. Surface Temperature vs. HCMM Temperature 4 5 Offset 
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data occurred on the I3th day of May, while the ground surface data had been 
collected on both the 9th and 12th, the delta day values being 4 and 1, 
respectively. 

An examination of various atmospheric parameters corresponding to 
this time period yielded some significant observations. The average air 
temperature for this period had risen 14.5 °F, while the discorded high pan 
temperature had risen 21 °F from the 8th to the 10th, dropped 6 °F on the 11th 
and risen another 7 °F by the 13th. Based on this data, it was detemlned 
that there had been a significant warming trend during this period, and that 
this trend would probably have been sufficient to have a marked Influence on 
the lake surface temperature. A rigorous examination of the data that was 
collected on days other than HCMM overflight days was made In order to Insure 
that circumstances such as this did not Influence the calibration of the HCMM 
equation. 


8.2.2 Detailed Ajialysls 

A plot was made of the yearly record of sever/il atmospheric 
parameters that were felt to be Indicators of conditions on the lake surface. 
Table 8.1 Is a list of the parameters used, as well as the sources of the 
data. Care was taken to select climatological data from stations felt to be 
representative of the study area, but In some cases the best data obtainable 
were still subject to question. The source of %he climatological data was 
taken Into account whenever use was made of It. 

Data deemed questionable for reasons dealing with changing 
atmr^spherlc conditions were not representative of conditions when the HCMM 
ovel'pass occurred, and therefore. Inappropriate for use In the calibration 
procedure. Such data were not totally discarded, however. Although not 
acceptable for the calibration procedure, this data was In some cases used In 
other phases, once the calibration had been completed. 

The determination of which periods of time were in need of attention 
was made. It had previously been determined that a plus or minus two degree 
variation In the ground data should probably be expected, due to the methods 
used In Its collection. In order to further refine this estimate, an 
examination was made of the variation of the residuals for the data taken on 
the same day as one of the HCMM overflights . A plot was made of the magnitude 
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Table 8.1 

Atmospheric Parameters Examined 


Parameter 

Source of Data 

Humidity 

Provo Airport, U.S. Bureau of 
Reclamation Weather Station 

Wind 

Utah Lake, Lehi NOAA Weather Station 

Sol-a-meter 

(sunlight intensity and duration) 

Provo Airport, U.S. Bureau of 
Reclamation Weather Station 

Evaporation 

Utah Lake, Lehi NOAA Weather Station 

Lake Stage* 

U.S. Bureau of Reclamation Record at 
BYU Boathouse 

Air Temperature 

Utah Lake, Lehi NOAA Weather Station 

High 


Low 


Average 


Water (evaporation pan) temperature 

Utah Lake, Lehi NOAA Weather Station 

High 


Low 


Average 


Cloud Cover 

(minutes of sunshine) 

Salt Lake City, Utah, National 
Weather Service Forecast Office 

*Lake Stage measurements were obtained by performing a linear Interpola- 
tion between. end of month readings. 
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of the residual for data collected on such days (delta day * 0). From this 
plot it was determined that the total spread is approximately 2.2 °C, somewhat 
smaller than what was initially expected. 

The average residual for each day of data collection was then 
calculated, and these residuals compared to the spread of same-day collected 
data. Any days with average residuals falling outside the spread of the same 
day data were then closely examined for possible reasons for the deviation. 

In most cases, atmospheric changes corresponding to the dates in question 
provided sufficient reason for the deletion of the data. It should be 
emphasized, however, that no data was deleted unless sufficient evidence was 
found to Justify the deletion. One notable example of this, occurred on the 
4th and 5th of September, 1979. This period had one of the largest deviations 
found, yet an extensive search of all records showed that it corresponded to a 
time of atmospheric stability. No rational whatsoever was found for its 
deletion. Consequently, it was left in the data base used for the calibration 
process, figure 8.11 gives the respective spreads in residuals for all data 
eventually deemed acceptable for use in the calibration procedure. 

8.2.3 SSR Results 

Figure 8.12 shows the SSR verses offset curve for all data not 
previously eliminated due to unstable atmospheric trends. As shown, the 
offset required based on this data is approximately +4.7 °C. This is in 
close agreement with the offset determined by several other researcher as 
noted in NASA (1980). 

Required offset vs. SSR curves were constructed for various subsets 
of the same data used in Figure 8.12. Figures 8.13 through 8.16 depict these 
curves for data obtained at the various time differences from the HGMM 
overflights (delta day ■ 0, +1, -I, and the combined -4 and -5, respectively). 
Table 8.2 summarizes all the results with September 4 and 5 data (discussed 
above) deleted. A summary of the results obtained for only same day data is 
found in Table 8.3. 

Upon examination and comparison of these two tables, it is evident 
that all required offsets are within the previously established spread of + 2 
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Figure 8.11. Total Spread Encountered in Average Residuals 
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Table 8.2 

Suimna ry of Required Offsets with Early 
September Data Deleted 


Subset 

Offset 

Number of Data Points 
Involved 

All Data 

4.89 

179 

North Lake 

5.16 

26 

Central Lake 

4.06 

24 

Provo Bay 

4.76 

48 

Goshen Bay 

5.13 

81 

July 

5.07 

30 

August 

5.05 

70 

September 

5.02 

49 

October 

4.06 

6 

November 

4.16 

24 

Same Day 

4.90 

61 

+ 1 Days 

4.80 

53 

- 1 Days 

4.95 

45 

- 4 Days 

5.47 

11 

- 5 Days 

4.47 

9 

Fish & Wildlife 

4.81 

137 

USBR 

5.17 

42 
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Table 8.3 

Summary of Offsets Using Only Same Day Data 


Subset 

Offset 

Number of Data Points 
Involved 

All Data 

4.90 

61 

North Lake 

4.61 

3 

Central Lake 

4.17 

8 

Provo Bay 

5.44 

19 

Goshen Bay 

4.77 

31 

July 

4.95 

8 

August 

5.50 

25 

September 

4.16 

12 

October 

4.04 

4 

November 

4.63 

12 



^C. In fact, assuming 4.7 to be the average offset and calculating to the 
nearest tenth of a degree, all data fall within +1 °C of the average. For 
this reason It was felt that the results obtained from this analysis were 
highly reliable, being based upon such a large amount of data with such a 
small relative variation. Table 8.4 summarizes the pvLvUn-'y tables Including 
the September data (all data) and also gives standard deviations. 

Based upon the above procedure and Information, It was determined 
that the proper offset to apply to the HCMM model Is 4.9 ‘’c. Although the SSR 
plot for all data used Indicates an offset of 4.7 It was felt that the 4th 
and 5th of September data was having a significant Impact, and although 
Insufficient evidence could be found to delete It from the data base. It 
differed from all other data sufficiently to warrant exclusion from the actual 
calibration procedure. 

8.3 Atmospheric Effects 

In order to evaluate the effect of various atmospheric parameters on 
the calibration, the value of the atmospheric factor being examined required 
to force agreement between surface measured and satellite measured 
temperatures was determined. Only same day data were analyzed. Table 8.5 
summarizes the correlation coefficients between the values of each Individual 
atmospheric parameter examined and the offset required. None of the 
parameters have a statistically significant effect on the HCMM calibration. 

However, the presence of a single point In the low air temperature 
range was all that prevented a high degree of correlation (0.740). Since this 
trend did not entirely disappear In either the average or high air temperature 
plots (both being considered better Indicators than the low air temperature). 
It was determined that a more extensive examination should be made. 

In an attempt to attain a better understanding of the air temperature 
effect all available data was analyzed, as opposed to only data collected the 
same day. The result of the addition of the supplementary data was to destroy 
the trends noticed when using the same day data only. There was not, however, 
a significant Increase in the amount of low temperature data. While It is 
possible that the addition of the "off day” data is masking trends. It appears 
that air temperature, as well as other parameters evaluated, does not have a 
measurable effect on the calibration of the HCMM equation over the range of 
temperatures examined. 
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Table 8.5 


Correlation Coefficients of Atmospheric Parameters 
and Required HCMM Offsets 


Parameter 

Correlation Coefficient 

Humidity 

-0.203 

Wind 

0.051 

. Sol -a -meter 

-0.247 

Evaporation 

0.074 

Lake Stage 

0.233 

Air Temperature 

High 

0.386 

Low 

0.486 

Average 

0.442 

Pan Temperature 

High 

0.240 

Low 

0.457 

Average 

0.339 

Cloud Cover 

-0.283 




8.4 Linear Model *' * » 

The development of a simple relationship between the available 
thermal Intensity and desired surface temperature values was considered. It 
was determined that a more accurate model (equation) could not be found, due 
to the theoretical backing of the HCMM model as opposed to the empirical 
nature of any other developed relationship. It was felt, however, that a 
computationally simpler model with similar accuracy could be developed. The 
relationship between temperature and emitted thermal Infrared radiation Is 
very nearly linear over small temperature ranges. A linear model for use over 
the anticipated range of surface temperatures could realize a significant 
savings In conversion and computation time. It should be emphasized that the 
linear model Is based on empirically determined constants and has no 
theoretical Justification. The model would only be useful If It were shown to 
be equally accurate as the HCMM supplied theoretical model. 

The procedure consisted of randomly dividing the data Into two 
groups, group A and group B. Group A was then used to re-callbrate the HCMM 
equation, while group B was used to determine the coefficients for the linear 
model. The re*-callbratlon of the HCMM equation was done In essentially the 
same manner as the Initial calibration, while the linear model was developed 
through the use of the Mlnitab statistical package which uses standard simple 
linear regression techniques. 

The re~callbrated HCMM equation was then applied to the data of group 
B and the SSR between actual surface and model predicted temperatures was 
determined. The linear model was then applied to group A data and the 
corresponding SSR value determined. 

With the SSR values determined through the above procedure, the ratio 
of the two values follow a statistical F~dlstrlbutlon, which allows the two 
models to be tested for equality of accuracy. The linear model developed 
using all available data is: 

T - -7.172 + 0.37494 I 

where T * temperature °C (0 < T < 30) 

and I “ HCMM thermal Infrared intensity 

Figure 8.17 shows the correlation between the linear model and the 
HCMM model in the range from 1 °C to 27 °C. The correlation coefficient of 
the above equation Is 0.982. It is suggested that the above equation can be 
used In cases where the utilization of a simpler model Is justified. 



Figure 8.17. Linear Model vs. HCMM Model +4.9 
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9. WATER QUALITY 

9*1 Methodology 

The one-to-one correlations presented in this section were developed 
with the HCMM data and water quality parameters according to the methods 
explained In Section 7. Correlation coefficients and graphs of the regression 
lines representing several data sets are Included. In each case the 
regression analysis shown was performed on a data set which represented the 
typical correlation of that particular comparison. A description of each of 
the parameters correlated In this analysis has been Included. 

9 .2 Day-Infrared Correlations 

The day Infrared (DIR) Intensity nearest each lake sampling station 
was Identified and tabulated for this correlation analysis. The net plankton 
measurements taken on Utah Lake have been used to represent the algae 
concentrations present In 1978. The major reason for this Is that more 
complete records were kept In 1978 of the net plankton counts on the lake than 
of the total or nanno plankton counts. Since the accuracy of this 
Investigation Increases with an Increase In the number of correlations made, 
the net plankton counts were chosen. These counts were assumed to be 
representative of the total plankton population In Utah Lake. They are 
hereafter referred to as plankton counts or plankton concentrations. 
Correlation of the total plankton counts was used for measurements made In 
1979. 

The correlation during the summer months of the log total plankton 
(algae) concentration In Utah Lake and the emitted heat (DIR) from the water 
surface Is very high for a naturally occurring system. Correlation 
coefficients are given In Table 9.1. The level of significance oc for each 
correlation Is also given. If the level of significance of the correlation 
between the algae concentrations In the lake and the DIR were 0.03 this would 
mean that all but 5 % of the distribution of the correlation could be 
Interpreted as being significant In the substantiation of this relationship. 

If the significance of the correlation were greater than 0.05 the correlation 
was considered "not significant" and was designated by n/s. 
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Table 9.1 

Correlation Coefficient of DIR and log Total Plankton 

1978 


Date 

0'7/28 

08/03 

08/10 

08/17 

08/23 

09/13 

10/14 

10/25 

Corr. 

0.341 

0.749 

0.640 

0.977 

0.783 

0.605 

0.758 

0,479 

Sig.(oc) 

0.005 

0.01 

0.025 

0.005 

0.005 

0.05 

0.005 

n/s 
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Table 9.2 

Correlation Coefficients of DIR and log of 
Ceratium Hirundinella Concentration 


Date 

07/28 

08/03 

08/10 

03/17 

03/23 

09/13 

Corr. 

0.754 

-0.224 

0.220 

0.533 

-0.176 

0.705 

Sig.(cv) 0.01 

n/s 

n/s 

n/s 

n/s 

0.025 
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As the algae concentration In the water increases It appears more 
like a land mass on the DIR printouts; l»e., collecting more he^s^t during the 
daytime hours so that It appears warmer and then losing ,heat more rapidly than 
the clearer waters during the nighttime. This phenomenon'. can be observed very 
easily by comparing the day and night Infrared data as d|.splayed by*j|^ c^lor 
graphics. The color pictures of the DIR and NIR (night Infrared) data clearly 
show the same areas In Utah Lake being respectively warmer and cooler. 

The description of the waters concentrated with algae being similar 
to land bodies Is valid since the same phenomenon Is responsible for the 
behavior of each with respect to temperature. The addition of the algae to a 
mass of water has the effect of lowering the specific heat of the mass from 
the specific heat of clear water (1.00 cal/gm-’^C). Since the specific heat of 
land Is also lower than 1.0, both the land and the algae concentrated water 
acquire and lose heat more rapidly than clear water. This makes possible the 
detection of algae concentrations In water bodies such as Utah Lake by . 
monitoring the heat from the lake during the day and night. Figure 9.1 Is a 
graphical representation of the linear regression of the DiR and log total 
plankton as measured on Utah Lake, July 28, 1978. The regression equation and 
the value of r are Included. 

Ceratlum hirundlnella (Figure 9.2) Is a solitary, unicellular plant. 
Its cells are generally narrowly fusiform with one apical horn and 2 or 3 
stouter and shorter basal horns. The cells are 30 to 72 ym wide and 100 to 
400 ym long. Although rare In early June on Utah Lake, this algae was one of 
the dominant plankton throughout the rest of the summer. It was often 
abundant enough to color the water muddy-brown and plug plankton nets 
(Rushforth, et al, 1980). Correlation and significance levels for CH are 
given In Table 9.2. 

No real, continuous relationship could be established between the 
Ceratlum hirundlnella population In Utah Lake and the emitted heat (DIR) from 
the water surface. The most probable reason for this Is that the Ceratlum do 
not form clumps or clusters which would Increase their density In the water 
and thereby Increase both heat conductivity and storage. Another Important 
observation Is that while Ceratlum is the second most dominant algae In the 
lake It constitutes less than 10^ of the total algae- count In the later summer 
months. This results In the Ceratlum not having much effect on the observable 
heat emlttance compared to the entire plankton mass. 







Aphanlzomenon flos-aquae (Figure 9<3) is a parallel trlchomes algae, 
united In bundles or flakes to form macroscop:^c agp;regates. Its apices are 
broadly rounded with cells 5 to 6 ym In diameter and 6 to 8 yro long. 
Aphanlzomenon contains numerous pseudovacuoles and oblong or cylindrical 
heterocysts. This blue-green algae was usually the most abundant and 
conspicuous summer plankton In Utah Lake. Table 9.3 shows correlations and 
significance levels for DIR and AFA. 

Figure 9*4 shows an obvious relationship exists between Utah Lake DIR 
In the summer and the Aphanlzomenon concentration In the lake. These 
correlations closely parallel those of total plankton which was expected since 
nearly 90% of the plankton count In Utah Lake In the late summer Is 
Aphanlzomenon. This study did not attempt to explain the reasons for the 

I 

change during the summer sea£:on from a large diversity of species of plankton 
to essentially only two by late summer. However, this phenomenon has been 
substantiated by the similarity of the Aphanf^yomenon and total plankton 
correlations with the DIR measurements. 

9.3 Reflectance Correlations 

The HCMM visual reflectivity measurements (W) were made by the 
satellite sensor and tabulated for this comparison. The VR Intensity nearest 
each station where lake measurements were made was recorded. Table 9.4 gives 
the correlations and significance of VR and total plankton relationships. 

During the early summer and on through the middle of August a very 
high negative correlation exists between the algae concentration In the lake 
waters and the reflectivity of the lake surface. This negative correlation 
Indicates that as the concentration of the algae Increases In the Utah Lake 
waters the reflectiveness of the water surface decreases or appears darker. 
Figure 9.5 Is a representation of the regression of the VR and log Total 
Plankton on August 17, 1978. 

In late August this Inverse relationship between VR and total 
plankton concentration decreases and by October there is essentially no 
relationship. The reason for the breakdown of this relationship appears to be 
the occurrence of blooms of algae on the lake surface. These whitish, 
floating matts of algae and air bubbles are very reflective. Development of 
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^ * Table 9.3 

Correlation Coefficients of DIR and log 
Aphanizomenon Flos-Aquae Concentration 


1978 


Date 

07/28 

08/03 

08/10 

08/17 

08/23 

09/13 

10/14 

10/25 

Corr. 

0.993 

0.728 

0.628 

0.987 

0.783 

0.599 

0.767 

0.503 

Sig.(flf) 0.005 

0.01 

0.025 

0.005 

0.005 

0.05 

0.005 

n/s 


Table 9.4 


Correlation Coefficients of VR and log 
Total Plankton Concentration 
1978 


Data 

07/23 

08/03 

08/10 

08/17 

08/;?3 

09/13 

10/14 

10/25 

Corr. 

-0.960 

-0.946 

-0.972 

-0.979 

-0.652 

-0.740 

-0.190 

0.005 

Sig.fcl'} 

0.005 

0.005 

0.005 

0.005 

0.025 

0.025 

n/s 

n/s 


DIR (HCMM INTENSITY UNITS) 










these reflective blooms over the area of highest algae concentration which had 
previously been the least reflective areas on the lake effectively reduces the 
Inverse correlation between the HCMM VR and the algae concentrations. 

While some correlations between Che VR and log of net Ceratlum 
hlrundlnella count are quite significant, no continuous relationship could be 
established. The probable reasons for this are Che same as those for which no 
continuous relationship could be established between Ceratlum and DIR; l.e., 
the Ceratlum cells do not Join together to form larger bodies In the water and 
the proportion of the total plankton concentration made up by Ceratlum Is 
relatively small. This lack of correlation and significance are shown In 
Table 9.5. 

Because the correlation analysis for VR and Aphanlzomenon shown In 
Table 9.6 doesn't give a different pattern of results than that given for VR 
and total plankton, a conclusion could not be reached about the reflective 
properties of Aphanlzomenon as compared with the total phytoplankton In Utah 
Lake. The decline In Che Inverse correlation beginning In late August Is 
caused by the same phenomenon, algal blooms, which Is responsible for the 
decline In the relationship between VR and the total plankton concentration. 
Figure 9.6 shows a graph of the linear regression of VR and the log of 
Aphanlzomenon flos-aquae concentrations In Utah Lake, August 10, 1978. 

9.4 Nutrients 

Analyses were made on samples of Utah Lake waters taken at the lake 
measuring stations to determine the amount of phosphoitts present In Che 
sample. These analyses were made in the Brigham Young University Water 
Research Laboratory. The measurements of phosphorus concentration which 
coincide with Che algae and HCMM measurements made on Utah Lake have been 
compiled by date of sampling. Table 9.7 gives Che correlations and 
significance of the phosphorus and plankton concentrations. 

These correlations Indicate that at least In the first half of the 
year the plankton concentration In Utah Lake Is dependent on the concentration 
of phosphorus present. This was to be expected since phosphorus Is a basic 
nutrient of nearly all blue-green algae. There were only two dates on which 
data was available to analyze the correlation between the phosphorus 
concentrations and Aphanlzomenon concentrations In Utah Lake. Therefore, It 
Is difficult Co determine any meaningful relationship from these results. 
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Table 9.5 

Correlation Coefficients of VR ano log 
of Ceratium Hirundinella 
1978 


Date 07/28 

08/03 

08/10 

08/17 

08/23 

09/13' 

Corr. -0.883 

-0.362 

-0.467 

-0.571 

-0.313 

-0.5£.= 

Sig.(a) 0.005 

n/s 

n/s 

n/s 

n/s 

0.05 


Table 9,6 

Correlation Coefficients of VR and 
log .Aphanizomenon Concentration 


1978 


Date 

07/28 

08/03 

08/10 

08/17 

08/23 

09/13 

10/14 

10/25 

Corr. 

-0.945 

-0.836 

-0.949 

-0.994 

-0.645 

-0.731 

-0.214 

-0,005 

S1g. (Of) 

0.005 

0,005 

0,005 

0.005 

0.025 

0.025 

n/s 

n/s 


ORIGINAL PAGE IS 
OF POOR QUA15TY 


Y-- S2.38 - 5.990X 


X= LOG OF AFHANfZOMENON FLOS-AQUAE CONCENTRATION 

(NUMBER /LITER) 

Figure 9.6, Plot of VR and log of Aphanizomenon flos-aquae Concentra- 
tions on Utah Lake, August 10, 1978. 












Table 9.7 

Correlation Coefficients of Phosphorus and log 
of the Total Plankton Concentrations 



DIR and 

loq Total 

Plankton 

DIR and loo 

Afanizomenon 

Date 

07/2S/78 

03/10/79 

05/12/79 

07/28/78 

05/12/79 

Corr , 

0.488 

0.566 

0.725 

0.864 

0.384 

Stg. 

n/s 

0.05 

0.01 

0.005 

n/s 


Table 9.3 

Correlation Coefficients of Nitrogen and the log 


ui inti 

Date 

1 uua 1 r 1 OMN 

07/23/73 

03/10/79 

a L 1 uno 

05/12/79 

Corr. 

0.688 

0.556 

-0,066 

Sig. fcl') 

0.025 

0.05 

n/s 


Analyses were performed in the Brigham Young University Water 
Research Laboratory on water samples taken at testing stations on Utah Lake in 
order to also determine the concentration of nitrogen In the samples* The 
nitrogen present was classified as ammonia, nitrite or nitrate* The total of 
the nitrogen measurements which coincide with algae and HCMM measurements made 
on Utah Lake have been compiled by date of sampling* Corr/ilatlons of these 
measurements are shown In Table 9*8. 

Nitrogen Is a nutrient necessary for algae growth; however, the 
results of this correlation analysis are Inconclusive. This Is because there 
were only a few days on which measurements were made of both nitrogen and 

algae concentrations. 

» 

9 .5 Turbidity 

Turbidity is the measure of opacity or light scattering properties of 
the water sample. Turbidity measurements of several Utah Lake water samples 
have been made with a Model 2100A Turbidimeter in the Brigham Young University 
Water Research Laboratory. This meter measures In Nephelometric Turbidity 
Units (NTU) which are a measure of the Intensity of the light scattered In the 
sampl at an angle of 90° from the direction of the incident light source. 

The readings which coincide with the algae and HCMM measurements made on Utah 
Lake have been compiled by date of sampling. Correlations of these 
measurements presented In Table 9.9. 

There does not appear to be any significant correlation between the 
turbidity of the water samples and the concentration of the plankton In them. 
This may be explained by several factors. First, Utah Lake Is a very shallow 
lake so winds cause complete turnover In the lake waters and stir up Into 
suspension the bottom sediments. These suspended sediments may be principally 
responsible for the turbidity measured In the water samples. Secondly, the 
algae present In the water column are not so likely to scatter Incident light 
as suspended solids because they are relatively transparent when not clustered 
together and allow more light to be passed or refracted than diffracted. 

Correlations from Table 9.10 show that the more turbid waters of Utah 
Lake do not appear to emit any more heat (DIR) than the other waters of the 
lake. This Is consistent with the result that there was no significant 
correlation between algae concentrations and turbidity. 
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Table 9.9 

Correlation Coefficients of Turbidity and the log 
_^^^ofth|TotaT_^_nankton^_Conc^ 


Date 

07/28/78 

03/10/79 

Corr. 

-0.121 

0.520 

Sig. (cy) 

. n/s 

n/s 


Table 9.10 

Correlation Coefficients of DIR 

and Turbidity 

Date 07/28/78 

07/10/79 

08/24/79 

09/22/79 

Corr, 0.158 

-0.150 

0.240 

-0.456 

Sig.ior) n/s 

n/s 

n/s 

n/s 


Table 9.11 

Correlation Coefficients of VP. and Turbidity 


Date 

07/28/73 

08/31/78 

07/10/79 

08/24/79 

0911217 ^ 

Corr. 

0.137 

-0.189 

0.604 

0.090 

0.857 

Siq. 

ri/s 

n/s 

0.025 

n/s 

n/s 
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As Table 9*11 shows, there does not appear to he any significant 
continuous relationship between the reflectivity of Utah Lake and the 
turbidity of the water. Again, this confirms the results of the correlation 
between plankton concentrations and turbidity; l.e., that there was no 
significant relationship between them. 

9.6 Day Infrared and Reflectance Correlations 

The relationships of Table 9.12 show a certain level of negative 
correlation which was expected in this analysis since the DIR had a strong 
positive correlation with algae concentrations while VR had an even stronger 
negative correlation with the same. These results are simply a statistical 
expression of the fact that algae may be presumed to be present in the waters 
of Utah Lake wherever VR intensities are relatively low. 

From early summer until late fall the DIR and VR color graphics xrere 
very useful in Identifying areas within the lake of higher algae 
concentrations. The high negative correlation between the HCMM VR data on 
Utah Lake and the plankton concentrations Justifies the use of the VR color 
prints to locate algae concentrations. Areas which appear to be warmer on the 
DIR color prints and darker on the VR prints are areas of higher algae 
concentrations in the lake. Either the DIR prints or the VR prints could be 
used individually to predict the location of plankton concentrations In the 
lake, but the use of the two measurements together Increases the probability 
that the areas on the lake where both indicate algae are in fact locations of 
higher algae concentrations. 

Figures 9.7 and 9.8 illustrate this application of DIR and VR color 
pictures of the HCMM data. The DIR (Figure 9.7) and VR (Figure 9.8) sensed 
from Utah Lake on July 6, 1978 indicate a particularly high concentration of 
plankton in Provo Bay and a significant concentration in the middle of the 
lake. Reference to the tabulated measurement of plankton concentrations for 
July 6, 1978 substantiates the prediction of higher concentrations of algae in 
the areas Indicated. The highest recorded concentrations for the entire lake 
on this date were In Provo Bay. The stations along the eastern side of the 
lake near the n&ln body also indicated higher concentrations. 
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Table 9.12 

Correlation Coefficient!; of htr and VR 

Date 

07/28/78 08/31/78 07/10/79 08/24/79 

09/22/79 

Corr, 

-0,837 -0.796 -0,628 -0,535 

-0,547 

Sig. (a) 

0,005 0.005 0,025 0,05 

0.05 



Table 9.13 


Correlation Coefficients of NIR and log of 
the Total Plankton Concentration 

Date . 

08/31/78 

09/13/78 

Corr. 

-0.918 

-0.727 

Sig.(o<) 

0.005 

0.025 
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9.7 


Night Infrared Correlations 

The HCMM Night Infrared (NIR) measurements were made by the satellite 
and tabulated for this comparison. The NIR Intensity nearest each station 
where lake measurements were made was recorded. Table 9.13 gives the 
correlations and significance level of the NIR and plankton data. 

Due to the fact that only a small amount of NIR data was available 
for this study, only two dates on which NIR and total plankton measurements 
were made could be used. Data frasi j:#! few other days were acquired but were 
unusable because of clouds which made the Intensity readings Inaccurate. The 
two correlations which were performed however Indicate a strong negative 
relationship between the emitted heat from Utah lake at night (NIR) and algae 
concentrations In the lake waters. This result agrees with the observation 
made with respect to the correlation of DIR and plankton concentrations, that 
waters with significant suspended algae mass behave In a similar manner to the 
surrounding land. They heat up more quickly than the clear waters In the 
daytime and cool more quickly at night. Figure 9.9 shows a plot of the linear 
regression of the NIR and the log of total plankton count of Utah Lake, August 
31, 1978. 

Locating areas on the day and night IR color graphics scenes which 
are respectively warmer and cooler Is a method of Identifying areas within the 
laV.i with a relatively high algae concentration. Figures '^.lO and 9.11, the 
color representations of the HCMM DIR and NIR Intensities measured on Utah 
Lake September 23, 1978, Illustrate this relationship. The lake's northwest 
corner Is the warmest area during the day (Figure 9.10) and one of the coolest 
during the night (Figure 9.11), Indicating the presence of higher 
concentrations of algae. 
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10 . 


CORRELATION IMAGES 


10.1 Background 

The day infrared and visible (reflectance) imagery both show the 
effects of algae blooms and other phenomena. Warmer areas of the lake show in 
the Infrared. The reflectance data show the turbid areas. The Images 
separately however » give only a limited overall picture of the lake from which 
conclusions could be made. Therefore, a composite or correlation image was 
formed for each date. 

The study of the images showed that the algae blooms have high 
Infrared and high reflectance values. Other combinations of infrared, and 
reflectance intensities were evaluated in order to identify different unique 
features in various areas and throughout seasonal variations in the lake. No 
statistical studies were undertaken to relate the ground truth to the 
correlation image beyond the identification of seasonal patterns. 

HCMM coverage for 1978 and 1979 was sufficient to indicate definite 
patterns in the lake. Each part of the lake was analyzed separately 
throughout the entire year. Conclusions are given as to the cause of the 
patterns in the lake. The names of the months are abbreviated to the first or 
the first two letters except May is MY, March is MR, July is JL and September 
is S or SE. The number that starts the file name is the day of the month. 

10.2 Goshen Bay 

i 

Early summer and middle summer correlation Images are from March 
through July, as shown in Figure 10.1. The southwest side of Goshen Bay shows 
a low Infrared and a high visible reading (L-H) (refer to Section 7.4, Table 
7.1). This persists into the end of July. The 22 July picture (22 JL.COR) 
shows part of the Bay Interspersed with black, which because of its position 
With respect to the other colors is low Infrared and medium day visual (L-M), 
and dark blue (L-L) . The persistence of this pattern for several months would 
seem to overrule algae or diatom blooms as the cause. It is surmised that 
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this is the result of suspended sediments or turbidity caused by the wind. 

The south area of Goshen Bay Is very shallow and requires little wind action 
to cause turbulence. 

The northern end of Goshen Bay shows much more variation. During 
March It Is L-L and L-M. In May the boundary of Goshen Bay and the Center 
Lake area Is M-H, which Is yellow. This seems to Indicate a mixing between 
the Center Lake (M-M) and the north of the Bay (L-H) to produce this M-H 
Intensity. The July Images Indicate a return to the L~L and L-M Intensities 
that are normally associated with this part of the lake. 

The late summer sequence In Figure 10.2 follows the same pattern as 
that- of spring and summer. From the southwest to the northeast the change Is 
from L-H to a mixture of L-M and L-L. The Image of 11 August shows a 
disturbance In the southern part of the Bay. The reason Is unknown and the 
lake returns to the normal pattern. The September, late summer, Images 
confirm that the usual pattern In Goshen Bay Is disrupted and re-established 
throughout the summer (Figure 10.3). The variation In range of the 11 
September picture causes the masking of detail In the Central Lake as well as 
Goshen Bay. 

The fall Images In Figure 10.4 show much more variation In Goshen 
Bay. In dividing the infrared and visual Images Into three Intensity ranges, 
some detail may be lost because the large variation In the summer temperature 
and reflectivity ranges mask sub^‘ le differences. This contention Is not 
substantiated by the spring pi rvUtiis however, which are of a cold lake and 
show no masking of data due to variation In the Image. The final blooms and 
the death of the algae In October (Figure 10.4), probably account for this 
variation In Intensities along with some effect from the increasingly uniform 
temperature of the lake. 

After the 9 October Image the familiar pattern In the Bay disappears 
and Is replaced by a patchwork that Is of medium and high Infrared with the 
entire range of visible reflectivity high to low. The breakup of the algae 
has been indicated as a probable cause. In the 14 October Image there Is a 
possibility that this Is caused In some areas by groundwater from thermal 
springs. The northwest side of Goshen Bay is known to have a large number of 
springs. The temperature and amount of Inflow may be so small that they do not 
effect the temperature of the water enough to be seen until late autumn and 
winter when the lake has a cold and fairly uniform temperature. 



Figure 10.2 Summer Infrared - Visible Correlation Images. 
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Figure 10.3 Late Summer Infrared - Visible Correlation Images. 
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10.3 


Provo Bay 

The summer sequence, Figure 10.2, shows Provo Bay Initially with the 
Infrared varying from middle to high and the visible from middle to low. By 
the end of July this has changed to H-L and this Is maintained through the 
middle of August. By the 23 of August the homogenlety of Provo Bay Is gone. 

It Is largely a patchwork of M-M, M-L, L-L. With some disturbances and with 
varying contribution of different combinations, this pattern holds through the 
middle of October with the addition of some areas of M-H. 

The 25 October Image Is primarily L-M and L-H, which carries through 
to the 14 November Image, Figure 10.4. Provo Bay Is shallow and surrounded by 
swamps and reed beds. Combined with a constant fresh water Inflow, this 
causes Provo Bay to have a different environment. The decrease In water 
elevation would add to the Increase In reflectance. 

10.4 North Lake 

The north lake Is commonly thought of as the area which has the 
greatest biological activity. Because of this the comparison of Infrared and 
visible was expected to change radically and often. An Inspection of the 
Images throughout the year contradicts this assumption. The lake appears to 
return or attempts, to return to an equilibrium condition. This equilibrium 
condition Is one of a low Infrared and the full range of reflectance, from L 
to H. The north lake area seems to be where the disturbances In the lake 
migrate and not where they begin. It has been established that these 
disturbances as seen In the satellite Imagery are caused by biological 
activity. This has been confirmed In the study of algae. 

The winds affect the floating algae, Aphanlzomenon, and cause It to 
concentrate and to move. The monthly average of the wind direction at the 
Salt Lake City Airport is from 160 to 190 degrees. Thus, any floating material 
Is be pushed Into the north end of the lake. Inspection of the images 
Indicates that the disturbances start In the central part of the lake and 
migrate to the north. Conditions In the central lake and the prevailing winds 
seem to be the answers to questions of the conditions In the North of the 
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Center Lake 


On the 24 March Image the center of the lake Is predominantly M-M 
with large areas of M-H and M**L. This pattern fills all of the central 
portion of the lake from the east side to the west side. Throughout the 
season this same pattern forms In the central lake, migrates north or south, 
and then dissipates and reforms again In the center of the lake. This happens 
about once every month during the summer and the fall (Figures 10.1”10.4) . 

The major Inflow to the lake occurs on the south east side of the central 
section. This flow then goes north to the outlet. Most of the nutrients 
added to the lake would enter along the east side of the center of the lake* 

10.6 Algae Bloom Cycle 

There Is one sequence of color graphics that Is of particular 
Interest. Figures 10.5, 10.6, 10.7 and 10.8 show an algae bloom cycle which 
occurred during September 1979. They outline the changes that take place In 
the life cycle of an algae bloom. On each figure the top picture Is day 
Infrared, the middle Is the day visible and the bottom Is the Infrared-visible 
correlation. In the north central section of the first Image for 4 September 
1979 (Figure 10.5), there Is a small patch which displays medium Infrared 
(green, top picture), high reflectance (red and black, middle picture), and 
M-H correlation (yellow, bottom picture). This Is surrounded by areas of M^M 
and M-L. This patch Is the first Indication that alge Is concentrating. 

The 11 September Image (Figure 10.6) shows the bloom In full growth. 
At, this time the processes Involved In the transformation are not know. Many 
other Images of the lake show areas of the same M-H relationship and yet did 
not develop further into major blooms. The central part of the algae bloom Is 
a region of H-H. North and south of the center are areas of M-H and to the 
sides of this, M-M. The east-west edges of the bloom are of M-L Intensities. 
From this It is presumed that the greatest concentration is in the red (H-H) 
region with a decrease In concentration to the yellow (M-H)*and green (M-M) 
regions . 

The 16 September Image (Figure 10.7) shows the lake returning to the 
normal pattern. The warm area of 11 September has dissipated and the central 
lake avea is now the warmest. The temperature difference of the north lake 
between the 11 September and 16 September Images Is about 2 ^C. On 16 
September the north end of the lake exhibits a reflectance that Is only half 
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Figure 10.5 Day Infrared, Day Visible and Correlation 
Images, 4 September 1979. 
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Figure 10.6 Day Infrared, Day Visible and Correlation 
Images, 11 September 1979 . 
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Figure 10.7 Day Infrared, Day Visible and Correlation 
Images, 16 September 1979. 
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Figure 10.8 Day Infrared, Day Visible Correlation 
Images, 21 September 1979. 
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as gredt as the 11 September Imaget The change between the 11 September and 
16 September Is marked and Indicates a dispersal of the bloom that is equally 
as rapid as the build-up from A September to 11 September. 

The correlation image for 16 September shows the center of the lake 
as pink, M-L. The east side of the north lake has a M-H and M-M. |;0Ct@tn, 
which could be the algae remaining from the bloom# The boundary &£ the center 
lake and Goshen Bay displays a mixing pattern from M-L in the center lake to 
L-L in Goshen Bay. The algae bloom has been dispersed Over the north lake and 
that from the center lake to Goshen Bay there is a drop in water temperature 
across a mixing zone. 

The 21 September image (Figure 10.8) shows a more normal pattern. 

The only anomaly in the 21 September image is that where the algae bloom once 
was is a large area of low infrared, medium visible, and therefore L-M 
correlation. It is possible that the dying or dispersing algae are of such 
low concentration that there is little heat produced from the photosynthetic 
effect and that the relatively high reflectance causes less heating by solar 
radiation, resulting in a condition of low temperature and middle reflectance. 

The degree of statistical correlation between the reflectance and 
plankton concentration data decreased sharply from a very strong negative 
relationship during the early summer to practically no relationship by 
October. The reason for this loss of a statistical relationship appears to be 
that algal blooms begin to form in the late summer and Increase in number and 
size through the latter part of October. These blooms are characterized by 
large floating masses of algae which appear whitish due to the large number of 
entrained air bubbles produced by the photosynthetic r<rocesses of the 
plankton. These light colored floating masses are very reflective and reduce 
the relationship between alj^ae concentration in the water and the accompanying 
lower reflectivity of the earlier season. 
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EVAPORATION 


11.1 Dalton Equation 

i The evaporation analysla first Involved application of the ao called 
Dalton type equation. The typical form la: 

where E *■ evaporation 

eg“ vapor pressure at the evaporating surface. 

vapor pressure In the atmosphere above, 
f (u) ■ function of the horizontal wind velocity of the form 

a + bu, where a and b are constants and u Is the wind- 
speed. 

This equation was used because of Its simplicity and the ability to use 
satellite data directly In the calculation of the saturation vapor pressure 
term In the equation. 

Lake temperatures were calculated with the HCMM equation (NASA, I960) 
and given the predetermined 4.9 degree centigrade offset. The corrected 
temperatures were then used to calculate the saturation vapor pressure of the 
surface of the lake. Maximum air temperatures from the Lehl station were used 
liX calculating the prevailing vapor pressure. 

11.2 Wlndspeed Function 

{'/1th the corresponding Lehl pan evaporation, the values for the 
Wlndspeed function f(u) ■ E/(eg-e^i) were calculated. These values were 
plotted against both the Lehl and Provo Airport wind data. Insignificant 
correlation existed In both cases. A comparison of these preliminary plots 
with those of Penman (1948) showed that the wlndspeeds of the HCMM days did 
not have as great a range as those of Penman. This Is probably because the 
HCMM Images are only for days that are clear and relatively calm, being In 
between storms. Therefore days with high wlndspeeds did not have 
corresponding HCMM data and were not studied In this analysis. 
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other ways were explored In which the HCMM data could be used In the 
evaluation of wlndspeed functions. The relationship between maximum pan 
temperature and maximum HCMM lake'^^tempe^Tature Is seen In Figure 11.1 This 
figure shows the result of two wlndspee^ functions, from Lehl pan data and 
satellite data, being plotted on the same graph. 

In the computation of the pan wlndspeed function the average air 
temperature and average humidity were used to calculate the dewpoint 
temperature. From this the vapor pressure of the air, eg, was calculated. 

The maximum pan temperature yielded Cg (maximum pan), the water surface vapor 
pressure. Maximum pan temperature was chosen because It corresponds to the 
HCMM day measurement which was' also a maxlmuni temperature. These maximum pan 
wlndspeed function values, f(u), were calculated and plotted against the 
wlndspeed u In Figure 11.1. 

The maximum HCMM wlndspeed function was calculated with the e^ as 
above and Lehl evaporation* The water surface vapor pressure, eg (HCMM day), 
was obtained from the use of HCMM lake temperature data for the fourteen HCMM 
days. The wlndspeed functions were then also plotted on Figure 11.1. 
Superimposing the two wlndspeed functions In this manner reveals the 
differences. 

The wlndspeed function values from the HCMM data are much higher than 
those form the Lehl pan data. This Is caused by the difference In temperature 
of the lake and the pan. For the days studied the pan temperature ranged from 
seven to fourteen degrees centigrade hotter than the lake. This causes Cg for 
the HCMM values to be smaller and therefore gives a larger f(u). Analysis of 
this difference would be useful In the study of pan coefficients, but Is 
beyrad the scope of this study and would require much more frequent satellite 
coverage of the lake. 

The same method of calculation used In Figure 11.1 was used In the 
plotting of Figure 11.2. Evaporation and e^ are the same as in Figure 11.1. 
Howev^^r, the average pan temperature was used to calculate the water surface 
vapor pressure, eg In the average pan wlndspeed function. 

The average HCMM wlndspeed function Is calculated as above except 
that the average of the HCMM day and night (DIR & NIR) lake temperature Is 
used In the computation of eg and hence f(u). Figure 11.2 Is a plot of the 
average Lehl wlndspeed function and the average HCMM wlndspeed function on the 
same graph. The superimposed values show that the average lake temperature Is 
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Figure 11. 1 . Windspeed Functions f(u) (in/mb) from ifaximum Pan 
Temperature (°C) and from HCMM Day Temperature (OC) versus 
Lehi Evaporation Pan Windspeed U (mpd) . 


€ 







C\ 


f I m) (IN/MB) 


O AVE. PAN TEMP., f(u) 
rs 0.617 

Y si, 184 4-0.03375 X 

• HCMM DAY/NIGHT AVE. TEMP., iS) 

rs 0.553 

Y s 1.728 * 0.02749 X 




X= WINDSPEED, U (MPD) 


figure 1 1 . 2. Windspeed Functions f(u) (in/mb) from Average 
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slightly lower than the average The average of the five 

day/nlght HCMM temperatures yields ■WlSdSJSJ lOnctlon that plots In the same 
region of the graph as the Lehl pan data. The study of evaporation and 
wlndspeed functions should be made with the use of average HCMM temperatures. 

Vapor Pressures 

An interesting relationship Is shown In Figure 11.3. Three plots are 
superimposed on the same graph. The first plot Is of eg (maximum pan) derived 
from the maximum pan temperature versus eg (average pan) derived from the 
average pan temperature. The second set of points plotted is of eg (maximum 
pan) versus eg (HCMM day/nlght) calculated from the HCMM day/nlght (DIR/NIR) 
average temperatures. The final plot Is of eg (maximum pan) versus eg (HCMM 
day) from the HCMM' day Infrared (maximum) temperatures. 

A comparison of the three sets shows that the eg values for the HCMM 
average temperatures are lower than those calculated from average pan 
temperature. The vapor pressure Calculated from the HCMM day temperature 
corresponds very well with the vapor pressure calculated from the average pan 
temperature. It Is possible then to use HCMM day data Instead of the 
day/nlght average. 

The agreement between the average pan and the HCMM day values Is very 
useful. Other researchers h.ave shown that a pan 12 feet in diameter has the 
same evaporation as does a lake. This allows the use of the relatively 
abundant day Infrared satellite data for the study of evaporation. 

11 .4 Evaporation Model 

Miller (1980) gives comparisons of evaporation from several different 
pans. The pan coefficients were derived using the water budget for the years 
1970-1973. The pan evaporation was compared to evaporation calculated from a 
model developed by F.I. Morton. These data were plotted with respect to the 
HCMM Intensities. Model and pan evaporation was by month; therefore, monthly 
HCMM averages were calculated using the average lake Intensity from the HCMM 
Images. This data, for the years 1978 and 1979, was then plotted against the 
DIR and the NIR HCMM Intensity as seen In Figure 11.4 and Figure 11.5. 

The correlation for both day and night Intensities was excellent. 

The plot of monthly model and pan evaporation versus DIR showed that 93 and 94 
percent respectively of the variation In evaporation could be accounted for by 
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HCMM day Intensity. DIR versus daily evaporation gave correlations of only 
0.8 to 0.86 (see Table 12.1, Section 12) as compared to 0.97 for monthly 
values. As with other methods of evaporation estimation, the shorter the time 
Interval the less accurate are the results. 

The plot of evaporation versus NIR shows slmlllar results (Figure 
11.5). The HCMM night Intensity accounts for 95 and 96 percent respectively 
of the variation in model and pan evaporation. However, the increase in 
accuracy by using the larger monthly period for night infrared la very small, 
from r ■ 0.95 to r « 0.98. This small change in accuracy indicates that 
estimates of evaporatatlon for any time period much shorter than a month 
should be made with the use of HCMM NIR* For periods longer than a month the 
day intensity is as accurate. 
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12. CLIH/VTOLOGICAL ANALYSIS 

12.1 Parameters 

The HCMM data were compared to the climatological data of 1978 and 
1979 In order to analyze and document the relationship to the hydrology of the 
study area. HCMM temperature data Included day, night, day/nlght average, and 
day/nlght difference. Several parameters were evaluated to determine If any 
significant correlations existed. These Included HCMM day and night data to 
pan evaporation and to mlnlumum and maximum pan temperatures, and wind to 
evaporation and surface temperature, and surface temperatures In different 
parts of the lake to each other. 

12 .2 Pan Evaporation 

The four types of HCMM temperatures, day, night, day/nlght difference 
and day/nlght average, were plotted against the Lehl same day evaporation and 
the average Lehl evaporation for the two days before the satellite overpass. 

A typical plot Is shown In Figure 12.1. This Is of HCMM night temperature 
versus Lehl average evaporation. The correlations and regressions of all the 
evaporation relationships are given In Table 12.1. Heat Is a main cause of 
evaporation of water and the expected high correlation was observed. 

The data for Provo Airport evaporation In Table 12.1 shows a 
difference In correlation from that of the Lehl evaporation data. One reason 
Is that Provo airport evaporation Is an average evaporation. Measurements 
were taken only when the pan was filled, which varied from two to six days. 

The effect of this Is unknown and the Provo Airport evaporation data was 
largely discounted. 

The two day average of Lehl evaporation versus HCMM temperatures gave 
better correlations than the same day Lehl evaporation versus HCMM 
temperatures except In the case of the HCMM day/nlght temperature difference. 
The use of a pan coefficient does not have an appreciable effect on the 
correlation coefficient and has only a small effect on the regression 
equation. Two of the cases In Table 12.1 were of particular Interest In this 
study. They were HCMM temperature difference versus evaporation and HCMM 
night temperature versus evaporation. 

. . V 
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Xs LEHI AVERAGE EVAPORATION (IN) 


Figure 12.1, HCMM Night Temperatures (“C) versus Lehi 
Average Evaporation (in). 









Table 12.1. Correlacion Coefficient 
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The comparison of HCMft' teroperatu\:e difference versus evaporation 
shoved high correlation for the same-day Lehl evaporation hut a much lower 
correlation with the two day Lehl average* This is explained by the fact that 
the two day Lehl average evaporation masks the effect on evaporation of sudden 
temperature changes. The same day evaporation is more sensitive to 
temperature changes. HCMM data for this study were from August and September 
of 1978 and 1979. For these months there is a large variation in day and 
night temperature. 

The heat storage of the lake maintained the water temperature as the 
air temperature fell during September and evaporation decreased by more than 
half during the month. The high correlation of the two is expected but is 
probably a measure of the increasing temperature difference occurlng at the 
same time as the evaporation decreased in early autumn. 

The best results were obtained from HCMM night temperature versus 
average Lehl evaporation. The correlation coefficient was 0.95 (Table 12.1). 
With proper calibration of the water body it would be possible to generate 
good evaporation estimates from the night infrared temperatures. Night 
temperatures on the lake are near minimum temperatures. There are two 
possible reasons for the good results. First, the average Lehl evaporation 
effectively masks any sudden changes In pan evaporation due to the changing 
local weather. Second, the night temperature readings were unaffected by the 
surface heat gain that is registered by the day readings. This means that the 
night temperatures are a much better gage of the average temperature and 
therefore heat storage of the lake. If nijht Infrared data become common it 
could very well be possible to use those measurements for energy budget 
calculations. 

12 .3 Wind and Surface Temperature 

The surface of the lake Is where energy is absorbed and released and 
much of this energy Is released through evaporation. The air-water Interface 
Is acted upon by the wind and the heat flux caused by sunlight. With the 
satellite imagery the effect of wind on surface temperature and hence 
evaporation was investigated. Table 12.2 gives the results of the study of 
various HCMM temperatures versus the wind* Wind measurements were from a 
continuous recorder at a height of 11.5 m that was maintained by the Bureau of 
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Table 12. 2, Correlation Coefficients, r, and 
Regression Coefficients, a and b, for Wind, 

Evaporation and Temperature Relationships where 
Evaporation (in) ■ a + b x Wind (mph)» 24 Hour Average. 


Lehi Average HCMM HCMM HCMM HCMM 

Evaporation Day Night Average Difference 



r - .276 

r 

- -.002 

r » 

.533 

r •» 

.355 

r * 

-.617 

Wind (mph) 

a - 5.243 

a 

6.684 

a * 

6.319 

a » 

10.78 

a » 

24 .52 

24 Hr. Avg. 

b - 4.965 

b 

- -.00071 

b * 

1 .504 

b « 

1.074 

b « 

-2.308 


Table 12. 3, Correlation Coefficients, r, and Regression 
Coefficients, a and b, for Utah Lake HC>C-1 Temperature 
versus Air Temperature, where HCMM Day Temperature 
(Oq) • a + b Air Temperature (®C). 


Ta 

Max 

Ta Average 

Ta min 

Ta 

rain 

r » 

.868 

r - .818 

r - .693 

r » 

.886 

HG>Ci DAY a * 

1.215 

a « 4.046 

a - -9.037 

a » 

9.226 HCMM WIGHT 

b 

1.405 

b - 1.221 

b ■ 1.036 

b * 

.5751 
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RtclaniAtlont HourXy wind apeeda gtaph and averaged. for the 

24 houri bafore the HGMM image, 2 pCmM^A 2%^ typical graph, wind veraua 
HCMM day/nlght temperature difference, is ahown in Figure 12. 2, 

The correlation of wind to evaporation la amall, but algnif leant, and 
confirms what haa already been auggeated. That ia, wind effacta evaporation 
only allghtly on large laHea* The comparlaon of wind to HCMM temperaturea 
gave mixed raaulta. There ia virtually no correlation between wind and HCMM 
day temperaturea. One reaaon for thla la that the IICMM data ia only for 
perioda of relative calm before and after atorma. kt theae tlmaa the 
windapeed la uaually within the range of only four to aeven tnllea per hour* 

The night temperature ahowed a algnif leant correlation. However, the night 
temperature of the lake ia also very constant during July, August, and 
September. The combination of a conatant wind and temperature would give a 
aignlflcant correlation only true for the specific condltlona and not 
generally applicable to the lake. 

The effect of wind on changing temperatures was investigated using a 
plot of 24 hour wind versus the HCMM temperature difference between day and 
night. This is ahown in Figure 12.2. The high negative correlation of "0.617 
means that the regreaalon line haa a negative slope. That ia, there is a 
aroaller tamperaturo difference with increasing wind. Thirty-elght percent of 
the variation in the HCMM temperature difference is caused by wind. During a 
windy day the daytime temperature of the lake rises less than on a comparable 
calm day. The wind increases the amount of energy removed by advectlon. 
Radiation and mixing of the warm upper layer in the water column accounts for 
the remainder of the loss. Further studies could be made using several night 
images. The correlation of the thermal Images with extensive wind 
measurements' made simultaneously on and around the lake could be the basis for 
additional studies. 

12.4 Air Tamperature 

Tables 12.3, 12.4, and 12.5 show the results from comparisons of air 
temperature to HCMM temperature. In Table 12.3 the HCM>t day and night 
temperature is compared with the maximum, minimum and average overland air 
tempar»b'‘’f®s» fbe maximum air temperature correlates very well, as the HCMM 
day temperature was designed to be a maximum at the earth's surface. 
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Table 12. 4. 65S*i;eaftt;i9,n“CQemc r, ^nd Regression 

Coefficients, a and b, for Air Temperatures versus HCMM Day 
Temperatures for Sectors of the Lake, where Air Temperature 
(°C) ■ a + b X HCMM Day Temperature (°C). 


'4 


HCMM Day Temperature 


Utah North Center Goshen Provo 

Lake Lake Lake Bay Bay 



r 


.868 

r 

m 

.746 

r 

Mb. 

M 

.857 

r 

m 

.888 

r 

m 

.813 

Ta max 

a 

H 

1^215 

a 

m 

6.026 

a 


2.060 

a 

m 

1.372 

a 

m 

2.407 


b 

• 

1.405 

b 

m 

1.150 

b 

■1 

1.348 

b 

■ 

1*446 

b 

m 

1.249 


r 


.693 

r 


.563 

r 

* 

.682 

r 

m 

.756 

r 

m 

.703 

Ta min 

a 


-9.307 

a 

m 

-4.734 

a 

m'- 

-8.597 

a 

m 

-8.923 

a 

m 

-10.14 


b 

■ 

1.036 

b 

m 

.8009 

b 


.9902 

b 

m 

1 .038 

b 

m 

.9971 


Table 12.5. Correlation Coefficient, r. Regression Coefficients, a 
and b, and Standard Deviation,^, for Utah Lake Average HCMM Day 
Temperature versus HCMM Day Temperatures for Sectors of the 
Lake, where UL HCMM Average Temperature (°C) " a + b x Lake ' 

Sector Temperature (°C). 


HCMM Day Temperature 



North 

Lake 

Center 

Lake 

Goshen 

Bay 

Provo 

Bay 

Utah Lake 

r - .945 

r - .991 

r » .994 

r - .938 

Average HCMM 

a ■ 1.696 

, a - .5209 

a - .02350 

a = 0.8167 

Day Temperature 

b - 0.9001 

b - .9631 

b -1.045 

b - 0.8902 

Std. Deviation 

1.03 

0*^ .421 

<r- .353 

1.09 






Likewise, the minimum air temperature correlates well with the HCMM night 
temperature* Had air temperatures b|een taken over the lake even higher 
correlations would probably have resulted. 

Utah Lake was broken into Its main sectors and the HCMM temperatures 
of each sector were compared to the average HCMM temperatures of the lake. 
Table 12.4 shows the correlation of maximum and minimum overland air 
temperatures with the HCMM day temperatures of the sectors of the lake. The 
results displayed in Table 12.4 are an extension of that given In Table 12.3. 
The temperature differences among parts of the lake were consistent over the 
three months of HCMM data. This is also shown In Table 12.5. Figure 12.3 
shows the good temperature correlation between the center sector and the 
entire lake. The summer months that the HCMM data covered, July, August, and 
September, are the most productive months for Utah Lake. The close agreement 
In temperature among the different parts of the lake during the summer should 
hold during the rest of the year when the lake Is less active biologically. 
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13. ALGAE AND EVAPORATION 




* 

13.1 Hypothesis 

During the summer of 1981 research was conducted to study the effect 
of algae on evaporation. Class A evaporation pans were seeded with algae and 
evaporation, temperature and biomass measurements were taken. Results of this 
were used to determine the extent to which algae altered pan evaporation. The 
conclusions from this study might be applied to the lake to give an Indication 
of the effect of algae In the natural environment. 

There are several processes blocking or driving evaporation that were 
of Interest. The first consideration was that of the Increased pan 
temperature due to photosynthesis by the algae. This would elevate the water 
temperature. The reflectivity of the algae affects the amount of energy 
absorbed In the pan. The blocking effect of a large mass of floating algae 
was the third consideration. The photosynthetic effect and the Increased 
absorption were expected to increase evaporation noticeably from that of a 
clean class A evaporation pan while the blocking action of the algae was 
expected to be minimal. 


13.2 Location and Methodology 

The Bureau of Reclamation (USBR) evaporation pan Is located at the 
Provo Airport. The enclosure where the algae pans for this study were 
located, is adjacent to the USBR site. Therefore the wind, rain and 
evaporation data to be used In this study were available from the Bureau of 
Reclamation. 

The size of the enclosure allowed three pans to be placed Inside. 

The construction of the stands and their Installation were In accordance with 
U.S. Weather Bureau Instructions and matched that of the Bureau of Reclamation 
(USBR) pan nearby. By this method any error due to difference In Installation 
was eliminated and the USBR pan measurements was used as a control In the 
experiment. 


117 


Measurement of evaporation was made with standard stilling wells 
Installed In every pan and the use of a hook gage> the hook gage gave 
readings to the nearest thousandth of an Inch. The same hook gage was used 
for reading each pan, eliminating error In the use of different Instruments. 
The measurements were read and the pans were filled at 10 a.m. each day. 

Thermometers were used In two of the pans to record the high and low 
temperatures. The 10 a.m. check of evaporation was late enough In the day 
that the thermometers could be reset to a low enough temperature that the next 
nights minimum and the day's maximum would move the Indicators and give a 
correct reading. The thermometers rested In the bottom of the pan and were 
usually read in place and then lifted out and reset. When the algae was thick 
the thermometers had to be lifted out and read Immediately. The south side of 
the bottom of the pan was shaded during most of the day and was where the 
thermometers were placed In order to give on accurate reading of the water 
temperature. 

13.3 Sampling 

Evaporation from the pan occurs at the air-water Interface and the 
algae studied was a floating variety. The samples ware taken from the 
surface. Concentration of algae varied over the surface of the pan. In 
addition, the wind would on occasion cause the algae to concentrate on one 
side of the pan. Sampling was done with a one quart plastic bottle. The 
mouth of the bottle was lowered Into the water to a depth of about 5 mm and 
the bottle was drawn across the pan at a uniform rate. The pattern used to 
sample was that of several straight lines. This ensured that the sample would 
be representative regardless of the variation In concentration. On those days 
that the concentration of algae was uniform a triangular pattern was used. 

The samples were taken to the B.Y.U. Environmental Analysis 
Laboratory and frozen. They were later thawed and tested for biomass. Each 
sample was homogenized and a part put Into a vial and centrifuged. The 
supernatant was drawn off and the algae washed Into an evaporation dish. The 
water was driven out by heating to 103 for at least eight hours. The 
dishes were weighed and the biomass calculated as mg/1. 


13.4 Algae Growth 

The early summer algae growth was from an Innoculatlon. The pans 
were prepared by filling them with water and adding a nutrient solution. A 
standard Allen's nutrient solution was considered to be satisfactory (Bold, 
1978 and Starr, 1978). Allen's solution was added to the pan at the same time 
as the water and allowed to set overnight. A sample of Anabaena flos-aqua, 
another prominent algae species In Utah Lake, was Introduced Into one of the 
pans In order to obtain a unlalgal culture. The algae did not grow. A larger 
sample of Anabaena was reseeded Into the same pan while the other two pans 
were seeded with algae taken from Provo Bay. 

The object of this proceedure was to grow algae in the pan and take 
samples as It grew from a small to a large concentration. The water turned 
green and the algae mixed In the water column, but no large surface 
concentration was formed. The Anabaena flos-aqua grew even less than the 
algae In the pans seeded from lake water. During the later part of the summer 
when the lake was very active, algae was collected and the pans' were stocked 
with this lake algae once a week using a quantity large enough to cause a high 
surface concentration. 

13.5 Errors 

Several phases of the data collection could have Introduced errors 
Into the study. The sampling time of 10 a.m. was Just before the major 
evaporation of the day took place. Variation In the time that the evaporation 
in tits pens was sampled would give erroneous readings. However, due to 
careful monitoring the error resulting from the variation in the sample time 
was minimized. Any sample time that differed from 10 a.m. was noted and the 
evaporation of the USBR control pan was extended for that period from the 
contlnous forms. 

The use of the USBR pan as the control pan was In itself a major 
source of error. The pan was filled with water approximately once a week and 
the measurements of evaporation and rain were recorded on charts with a 
Continuous recorder. The varying water level of the pan effected the amount 
of evaporation. Calculating evaporation from the strip charts was fairly 
reliable except during periods of rain. Because of this, all data during 
times of precipitation was disregarded. 


Sampling and determination of biomass was standardized as much as 
possible t The collection of the sample was the least reliable proceedure. 

The mouth of the bottle had to be Immersed the same amount each time and the 
speed with which It was drawn across the water had to be held as constant as 
possible* To reduce the amount of error, the same size bottles were used to 
obtain all samples. Also the same person took all the samples so that 
variation caused by different techniques would not be present. 

The centrifuge used to concentrate the algae would sometimes leave 
small pieces floating In the supernatant* Care was taken In this process and 
finally another more powerful machine was used* Error Introduced In this 
proceedure was considered minimal* 

13*6 Results 

Temperature was known to be Increased by masses of algae In the lake* 
Pan temperatures were taken Infrequently In the USSR pan* A comparison of the 
temperatures of the pans showed that the pans with algae had a smaller change 
In temperature from minimum to maximum* The amount of evaporation during this 
period did not vary significantly from pan to pan. The lack of difference In 
total evaporation Indicated that the amount of algae In the pans had no 
significant effect on evaporation* The top of the pans were covered with 
algae when many of the samples were taken* Table 13*1 shows a histogram of 
the middle of the Interval of the log of the algae biomass concentration and 
the number of observations In each Interval* The range of biomass data Is 
large enough that any effect on evaporation would have been shown* 

The amount of evaporation did not depend on the amount of algae 
present* The separate effects of differing albedo, photosynthesis, and 
blocking could not be determined* Figures 13*1 and 13.2 show the change In 
evaporation versus the log of the algae biomass concentration using both the 
USSR pan and the B.Y.U. weather station pan on the B.Y.U. campus. The 
correlation coefficients were almost zero and the regression lines were 
virtually horizontal. 

The plots from the B.Y.U. and the USBR pans are similar except for 
the differing Intercepts of the Y-axes. From this It Is concluded that any of 
the effects of algae on evaporation cancel each other with no net effect on 
evaporation* The use of B.Y.U* data results In a different Intercept because 
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Figure 13. 1-. Difference in Pan Evaporation (in), BYU Weather 
Station Pan minus NASA Pan, versus the log of the Algae 
Biomass (mg/l) . 
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the BtY«Ut pen Is located ovetr" Smiles* from the lake* This shows |:hat though 

• « 

the location of the pan nay be different* the regression line Is virtually 
horizontal* showing no correlation of algae and evaporation In both cases* 
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14. THERMAL SPRINGS 

14.1 Ground Identification 

It Is estimated that approximately 20 percent of the floor of Utah 
Lake more than one kilometer offshore Is spiked with springs or seeps 
(Brlmhall, et al, 1976). A large portion of these springs, Including those In 
the Goshen Bay area, are either known or suspected to be thermal In nature. 
During the summer 1975 a seml-lntenslve study of Utah Lake was made by means 
of a sonarllke device. The acoustical reflection pattern obtained from 
numerous profiles Imaged the thickness distribution and character of the 
underlying sediments. Using these prof lies, It was possible to Infer the 
existence and location of underground seeps and springs (Dustin and Merritt, 
1980). Figure 14.1 depicts the location of suspected spring and seep ares? In 
the lake bottom. 

Two years later, In 1977, aerial thermal scans were made by the Idaho 
Air National Guard. An RF-46 aircraft was used, flying at 5,000 ft. (Dustin, 
1978). A large number of surface anomalies were present on the imagery, 
s'l'iowliig small areas of thermal discontinuity. Figure 14.2 shows the location 
of the suspected areas. 

An observation made during the HCMM study period that warrants 
mention was the marked deviation of two boat measured surface temperature from 
the dally average on the 10th of July, 1979. Generally, the data obtained 
from different parts of the lake followed a fairly predictable pattern, with 
Provo Bay being the wairmest area and Goshen Bay the coolest. On the above 
date, however, tiro stations In Goshen Bay were significantly warmer than the 
remaining part of the lake. This corresponded to an atmospherically stable 
period substantiated by the lack of deviation of all other points sampled that 
same day. Upon close examination. It was found that the area containing these 
points quite often corresponded to a dissimilarity In surface temperature as 
sensed by the HCMM. 

Figure 14.3 shows an enlarged view of the Goshen Bay area, also 
depicting the location of the two deviant points. The close proximity of 
these two points to the strongest suspected spring area In Goshen Bay raised 
the question of a possible connection between the two. For the springs to 
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Figure 14.3. Proximity of Deviant Data to Suspected Spring Areas >> 
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have this noticed effect on the HCMM Imagery, one explanation would be that 
they have a somewhat sporadic, surging nature, releasing moderate amounts of 
warmer water upon occasion and being relatively Inactive the remainder of the 
time. However, this was not the suspected nature of the spring, but rather it 
was felt to have an essentially continuous flow. 

Another possible explanation was the development of a thermal "plume" 
emanating from the spring area. In atmospherically stable and windless 
periods, this plume might be established and measurable. Periods of high 
wind, conversely, would tend to mix and distribute the warm water with that of 
the surrounding area, making It's detection from the surface Impossible. A 
review of existing wind conditions at that time showed that for several days 
preceding the lake surface temperature measurement, winds were relatively 
calm. The next few days, however, between the lake surface measurement and 
the HCMM overflight, showed a marked Increase In wind activity, which was 
sufficient to have a mixing effect on the lake surface waters. This explains 
why the warm spots didn’t show up In the HCMM data of 14 July 1979. 

14.2 HCMM Identification 

One of the most Interesting and useful applications of the color 
grap’.klcs of the DIR and NIR HCMM data has been the Identification of these 
thermal areas within Utah Lake. The late summer and fall NIR color Images and 
the late fall DIR Images were used In an attempt to locate areas within the 
lake where significant thermal sources exist. 

As the temperature of the lake water cooled in the late summer an 
area of consistently higher NIR intensities developed in the southwestern 
quadrant of the main lake body. Several other areas of two or three pixel 
groups consistently indicated warmer temperatures as well. However, the 
southwestern quadrant area Is by far the largest and most easily observed. 
Figure 14.4, the color representation of the HCMM NIR Intensities measured on 
Utah Lake for September 4, 1979, shows the location and intensity of this 
thermal area (green and yellow) . 

By November the daytime temperatures of the lake had dropped to a 
level where some thermal areas could be seen In the DIR color prints. By this 
time the algae In the lake have been removed by the cold evening temperatures 
or have been carried to the northern end of the lake by the wind or with the 
flow toward the Jordan River outlet. Observed warm areas (yellow and pink) 
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Figure 14.4* NIR on Utah Lake, September 4, 1979. 



Figure 14,5. DIR on Utah Lake, November 14, 1979. 
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are therefore most likely a result of sources within the lake. In Figure 
14.5, HCMM DIR for November 14, 1979, the warmer pixels are In the same 
general area which showed warmer waters In Figure 14.4. This apparently 
substantiates the presence of thermal sources. 


15. 


GROUNDWATER 


15.1 Data Availability 

Initially, this research was to be done using chiefly quantitative 
data, obtained from various Federal and private organizations having conducted 
studies within the study area. A more thorough analysis of the available 
data, however, led to the realization that the vast majority of previous 
Studies had been done on the three deep aquifers and had virtually Ignored the 
near 'surface water table. The only available data was from county percolation 
testing done In the application for building permits. This data. In addition 
to being somewhat sketchy, was also subject to some question as to accuracy. 

The distinct lack of ‘what was felt to be a significant volume of 
usable data prompted the decision to expand the originally planned research 
Into the qualitative area. The data search had turned up several well 
documented maps covering the study area and showing regions of 
characteristically high and low groundwater. 

15.2 Quantitative Groundwater Analysis 

The quantitative portion of the groundwater study was procedurally 
similar to the development of the linear model for predicting surface 
temperature. The groundwater model was not restricted to a linear fit, 
however, so more advanced regression techniques were needed* These techniques 
are still matrix operations, however, and were done with the aid of the 
Mlnltab IX statistical package. 

Since no base model for this type of system was found In the 
literature, a variety of different types of curves were tried, using the 
correlation coefficient as the chief Indicator of the best fit. After the 
best curve was determined, the groundwater data was also extensively 
examined , through the use of the residual plots. In order to locate any 
unaccounted for parameter affecting the model. 

The quantitative data was obtained from unpublished records of county 
percolation tests performed In the study area. Twenty-two tests were recorded 
of which seven were unsupervised by county personnel. During the study 
period, a procedural change was made that required a county official to 
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directly supervise the percolation testing* The reason for this change was 
the Inaccuracy of reported results when testing was not overseen by a county 
official. Although data obtained from unsupervised tests were not Immediately 
rejected, they were closely examined as to their Impact on any results 
obtained through their Incorporation into the data base* 

15*2.1 All Data 

Figure 15*1 depicts the relationship between depth to groundwater and 
the HCMM sensed surface temperature, using both supervised and unsupervised 
data* As can be seen, no obvious strong correlation exists. It appeared that 
a best fit curve would possibly be a third or fourth order polynomial. Only 
0*394 was the best correlation coefficient for the fourth order equation given 
below: 

d ■ a + bT** , 

Where d > depth to groundwater 
and T ■ HCMM temperatures 

It Is recommended, however, that polynomials of 3rd order and higher 
be avoided, and that models of a different type l.e. log, exponential, etc>, 
be used (Ryan, et al, 1976). After extensive analysis of several recommended 
models and combinations. It was determined that the best fit to the curve was 
given by the following equation form: 



The correlation coefficient of the above equation was 0*411* This 
model, however, has absolutely no theoretical basis* It simply provided the 
best fit curve to the available data. 

15*2*2 Separated Data 

The data base was next separated Into two sets corresponding to the 
supervised and the unsupervised test* Figures 15*2 and 15*3 show the 
respective relationships between depth to groundwater and the HCMM sensed 
surface temperature* Again, no strong correlation exists* 

By using the same approach as above, several types of curves were 
examined In an attempt to find the ”best fit". Several equations yielded a 
correlation coefficient of 0*37 for the supervised test data, the simplest of 
which was a straight line: 



Figure 15.1. Depth to Groundwater vs. HCMi'l Surface Temperature 
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Figure 15.2. Depth to Groundwater vs. HCMM Temperature for 
Supervised Tests 
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Figure 15.3. Depth to Groundwater vs. HCMM Temperature for 
Unsupervised Tests 


t. 



d ■ a + bT 


Surprisingly » a correlation coefficient of 0*678 Was obtained through 
the use of the following model with only the unsupervised test data: 

d - a + bT '»/3 

Possible explanations as "^o'why the more questionable unaupervlsed data 
yielded the superior correlation coefficient are discussed later* 

15.2*3 Seasonal Variation 

Non-vegetated areas are preferred for this type of study for a 
variety of reasons* Therefore, the mid-summer values were deleted froa. the 
groundwater data base, leaving only the spring (April and May) and fall 
(November) data* It was hoped that through this approach the majority of the 
vegetation season could be eliminated, or at least the effect minimized* 

Table 15*1 lists the off season groundwater data used In the above 
correlation* As can be seen, the data Is essentially from 3 HCMM overpasses* 
Of the nine measurements, four are from May of 1978, four from April of 1979 
and one from November, 1979* Although the last data Is quite late In the 
season» an examination of National Weather Service data for the area showed no 
Indication of snow on the ground* 

Figure 15.4 shows groundwater depth plotted against HCMH temperature 
for the non-summer data only. The linear correlation of the groundwater depth 
to HCMM sensed temperature Increased to 0*61, which for some natural systems 
1$ considered acceptable 

Groundwater depth was then correlated with the HCMM intensity value 
as opposed to the computed surface temperature* With an exponential curve 
applied to the data a correlation coefficient of 0.64 was obtained. This 
higher correlation which was possibly due to the slight curvature of the 
relationship between HCMM Intensity and HCMM Temperature, yielded the 
following relationship: 

d » 10.827e ° I (HCMM) 

Of the off season data, all depths of greater than 60 inches were 
measured in May, while those depths less than 60 Inches were the April and 
November values. Since ground temperature would naturally be expected to rise 
as the spring season progressed, a measurable rise In sensed surface 
temperature would be expected regardless of groundwater fluctuation. This 
trend would likely hav^e a significant Impact on the relationship to be 
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Figure 15.4. Depth to Groundwater vs. HCMM Temperature for 
Off Season Data 





dav«l0p«d. thorofore, two indopondont appcoacheo mtn taken In an attempt to 
remove Ita effect* Firsts on attempt waa made to factor the eeaeonal 
temperature variation out of the data» reducing It to atandardlaed, and 
therefore comparable | values* Second, the different overflight dates were 
examined Independently. 

15*2.3*1 Temperature Standardization 

The first step In attempting to factor out seasonal variation was to 
find a suitable parameter that' could be used to reduce the data to seasonally 
Independent values* The temperature of dry soil at each of the data sites 
would have been Ideal, but no such data existed* The most reasonable 
approximation to this data that could be located was the recorded 8 Inch soil 
depth temperature at the Salt Lake City Airport Weather Station* As this was 
sufficiently far away (SO miles) from the study area to be quite questionable, 
the average dally air temperature for the study area was also recorded* Both 
the soil and air temperatures were each subtracted from the HCMM predicted 
surface temperatures* The resulting values, reflected the soil surface 
temperatures while minimizing the possible Influence of seasonal variation* 
This procedure had the effect of generally moving the four higher 
points to the left In relationship to the lower values as shown In Figures 
15*5 and 15*6* Although any significant correlation was essentially 
eliminated. It Is felt that In subtracting the seasonal variation, the true 
ground water/surface temperature relationship was more accurately established* 

15*2*3*2 Independent Overflights 

The separation of data Into Individual overflight days totally 
eliminated the aeasonal variation* However, It reduced the number of data 
points for each set to four, making significant results difficult to obtain* 
The two groups of data are shown In Figures 15*7 and 15*8* 

Figure 15.7 shows an Slight Indication of higher surface temperature 
with deeper groundwater, but the spread of the data Is quite broad and results 
are somewhat questionable* Figure 15.8 shows much better correlation between 
deeper groundwater and higher surface temperature* A curve through the points 
further substantiates that the relationship between groundwater depths and 
surface temperatures may follow an exponential curve of form: 
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Figure 15.6. Depth to Groundwater vs, HCMM Temperature 
Soil Temperature for Off Season Data 
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Figure 15.7. Depth to Groundwater vs. HCMM Temperature 
for May, 1978 Data 
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Figure 15.8. Depth to Groundwater vs. HCfiM Temperature 
for April , 1979 Data 



15.3 Qualltaftlve Groundwater Analysis 
15.3.1 Methodology 

The qualitative portion of the groundwater study took a much broader 
approach. Figure 15.9, 15.10 and 15.11 depict sketches of the maps obtained 
from the Utah State Geological and Mineral Survey CvflcjvMS ) , Soil Conservation 
Service (SCS) Septic Tank Study, and SCS high water table study, respectively. 

The method employed consisted of examining several transects through 
the study area. Transect location was based on several factors. Emphasis was 
given to areas of rapid change In groundwater depth. Nonvegetated areas were 
preferred, but as the study site was chiefly agricultural, this requirement 
was difficult to fulfill* Special attention was given to non-homogeneous 
areas, well defined shallow water table areas surrounded by deeper water 
table, etc. One problem encountered was disagreement among the three maps as 
to the location of some of the critical areas. Generally, this was resolved 
by taking transects based on a specific map, as opposed to constructing an 
aggregate map or using multipurpose locations. In alii, ten transects were 
chosen. Transect locations on the groundwater maps are shown In Figure 15,12, 
15.13 and 15.14. 

After having determined the key locations, the chosen transects were 
plotted on the overlays of the study area. Due to the expanded east west 
scale of the HCMM line printer maps, distances along the transects were rather 
severely dlstroted, depending on their relative position on the overlay. An 
east vest transect for instance, appeared roughly three times as long as a 
north south transect of similar actual length. 

Also due to the skew angle difference between day and night Images, 
the day/night lengths of the same transect varied. This variation in length 
became obvious as the overlays were used to obtain the HCMM sensed Intensities 
from the line printer maps. As all data from the same transect was to be 
compared to the various groundwater map derived groundwater patterns, It was 
necessary to adjust the lengths of the varying transects to correspond to a 
specific base length. The longest transect was chosen as the base length. In 
using the longest. It was possible to Insure that no critical data would be 
deleted. As transects needed to be lengthened, the number of additional 
Intensity values were determined. The transects were then proportioned Into 





Figure 15.9. Utah Valley Groundwater Depth Contours According to Utah State Geological and Mineral Survey 







Figure 15.10. Utah Valley Areas of Seasonal GrounAKater Depth Less than Four Feet (U.S. Soil Conservation 
Service f-lap) 
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Figure 15. 11. Utah Valley Areas of Poorly Drained Soils (Utah County Government Map) 
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Figure 15 . 12 A. Transects Examined as Compared to USGMS GroundMater Map 




Figure 15.I2B. Transects Examined as Compared to USGMS GrounA<ater Map 
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Figure 15.13A. Trdnsects Exdinined ds Cocnpdred to USSCS Groundn^dter Mdp 
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Utah County Government GrouncWater Map 




Figure 15.14B. Transects Examined as Compared to Utah County Government Groundwater Map. 







that number of segments and additional Intensity values were Inserted between 
appropriate segments* The magnitude of the Inserted value was the arithmetic 
mean of the two neighboring values* 

After the ten best transects had been determined they were plotted as 
block diagrams, with the groundwater condition being the ordinate and distance 
along the transect being the abslclssa* These transects were then overlayed 
on the HCMM DIR Intensity maps for various seasonally representative days* The 
HCMM values were plotted as the ordinate, again with distance along the 
transect as the abscissa* These plots were then overlayed onto the block 
diagram of the transects previously obtained* 

15*3*2 Results 

A detailed comparison of the transects failed to lead to any 
conclusive results. Although some transects showed strong evidence of 
correlation, there existed a disturbing sign descrepancy. Le Shack et al 
(1975) had also experienced a similar discrepancy while doing work with soil 
moisture. Figures 15.15 and 15.16 represent two of the best correlations 
obtained* It was observed that as the groundwater approached the surface, the 
HCMM DIR Intensities tended to rise, Indicating higher ground surface 
temperatures* This result was opposite of what was expected, as well as 
opposite to the results obtained in the quantitative portion of the study* 
Apparently, what was being measured in the HCMM DIR Intensities was some 
effect other than groundwater* 

An examination of only early spring or late fall data failed to 
Improve the situation* Figure 15*17 shows the correlation of transect 1 for 
May of 1978* Although a correlation seems to exist. It appears to change sign 
from one end of the transect to the other* The reason for this was 
unresolved* 




I DEPTH TO GROUNDWATER (ft) 

Figure 15.15. Comparison of Groundwater Profile to HCMM Intensities 
for 14, July 1979 (DAY INFRARED) Transect #5) 
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16. 


SUMMARY 


16.1 Calibration 

The study results Indicate that the appropriate offset to be applied 
to calibrate the HCMM model to the Utah Valley area Is + 5 degrees Celsius. 

The use of this offset yields radiative surface temperatures within 
approximately 1/2 degree Celsius of ground truth. This variation corresponds 
roughly to the accuracy of the combined sensing and telemetry equipment. 

The lack of correlation between the required offset and the various 
examined atmospheric parameters Indicates that the + 5 degrees should apply to 
any region of similar climate. The only significant parameter not examined 
was atmospheric translmlsslvlty. This should be based on suspended 
particulate matter, and/or total atmosphere thickness above the study area. 

The 5 degree offset, therefore, may be restricted to areas of corresponding 
elevation (approximately 4,500 feet above mean sea level) and air quality. 

In cases where a simplified relationship between HCMM sensed thermal 
IR intensities and surface temperatures is justified, it is felt that the 
linear relationship below yields sufficiently accurate values over the 
specified temperature range. 

T - -7.172 + 0.3749 I 

where T « temperature in degrees Celsius (0 < T < 30) 
and I ■ HCMM sensed thermal I R value 

Although this is a fairly narrow temperature range, it essentially spans the 
range of water surface temperatures encountered In the study area. It should 
be eraphaslaed that this model is entirely empirical, having no theoretical 
basis. 

16.2 Water Quality 

The st«s\tlstlcal analyses performed on the data obtained from the Heat 
Capacity Mapping Mission and the water quality measurements made in Utah lake 
in 1978 and 1979 yielded several significant relationships. The strong 
correlations identified between the DIR, VR and plankton concentration 
measurements provided the foundation for the use of the color prints to locate 
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areas of algae growth In Utah Lake> Several methods for identifying algae 
concentrations In the lake were determined, Including the comparison of DIR 
and VR prints and the comparison of DIR and NIR prints. 

Correlations between turbidity and the HCMM data were Insignificant. 
This demonstrated that suspended algae In the lake were much more responsible 
for changes In the water's reflectivity than suspended sediments. 

Relationships developed between algae concentrations In the lake and 
nutrient levels were inconclusive. This precluded the possibility of 
estimating nutrient concentrations In Utah Lake by assuming a relationship 
with algae growth and HCMM data. 

16.3 Evaporation 

Comparisons of wind and wind speed functions derived from HCMM data 
failed to yield significant correlations. This was mainly due to the small 
range of wind speeds on fair weather HCMM overflight dates. The HCMM 
day/night average temperatures, and likewise wind speed functions, compare 
well with average evaporation pan temperatures and extrapolation from pan data 
for larger wind speed ranges can be made. HCMM day temperatures are much 
lower (7°to 14°) than maximum pan temperatures yielding much higher windspeed 
functions. Vapor pressures calculated from the HCMM day temperatures 
correspond very well with the vapor pressure calculated from the average pan 
temperatures. 

Utah Lake monthly modeled evaporation and pan evaporation correlated 
well with HCMM day and night temperatures. The correlation of daily 
evaporation with HCMM day temperatures was lower (0.8) than with HCMM night 
temperature (0.97). For short time intervals (days), the night data might be 
useful in estimating evaporation while for monthly evaporation either night or 
day data appears useful. The encouraging night data results are most likely 
because night temperatures are unaffected by surface heat gain that Is 
significant during the day. That Is, night temperatures are much better gages 
of the actual average temperature and heat storage of the lake. 


16.4 Groundwater 

The study relating off season groundwater depth to HCMM sensed 
thermal IR Intensity yielded a model with an extremely good correlation to 
ground truth. However, further research is needed before it can be accepted 
as significant. The limited amount of data is the problem. Also, the 
influence of the extraneous parameters mentioned in the analysis of the 
quantitative groundwater data could be of equal or greater magnitude than the 
actual data being studied. It is felt, however, that there does exist a 
measurable relationship between HCMM sensed surface temperatures and depth to 
groundwater. This relationship appears to follow an exponential curve of the 
form} 

d “ ae^^ 

where d “ depth to groundwater 

T “ HCMM sensed surface tevoperatures 
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